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PREFACE 


The  following  Technical  Report  Is  an  adaptation  of  an  M.S.  thesis 
submitted  to  The  University  of  Tennessee  in  August  1982.  It  Is  the 
result  of  an  effort  to  go  beyond  the  usual  Interpretation  of  a 
diffraction  field  and  to  Interpret  simultaneously  the  effects  of  linear 
diffraction  and  nonparallelism  In  terms  of  pulse  echo  amplitudes 
observed  on  an  oscilloscope.  It  is  a  necessary  first  step  toward  a 
more  general  interpretation  which  Includes  the  nonlinearity  of  the 
solid  propagating  medium. 

Even  in  the  linear  approximation  the  evaluation  of  the  amplitudes, 
and  hence  the  attenuation,  of  pulsed  ultrasonic  waves  in  solids  can  be 
complicated  by  a  number  of  factors.  Of  those  factors  most  often 
considered,  this  report  concentrates  on  extending  our  understanding  of 
diffraction  and  nonparallelism.  A  mathematical  model  for  correcting 
the  echo  height  for  the  effect  of  a  wedge-shaped  sample  is  constructed. 

By  using  a  unique  coordinate  system,  the  path  of  a  multiply-reflected 
ultrasonic  wave  is  transformed  into  an  equivalent  unidirectional  path. 

A  dlffractton  correction  then  is  applied  by  numerically  Integrating 
an  improved  version  of  the  farfleld  solution  to  the  diffraction 
Integral.  This  model  Is  used  to  interpret  data  taken  on  a  steel  plate 
which  has  different  facets  ground  over  a  range  of  angles  between  0  and 
0.01175  radians.  Plots  are  given  of  corrected  attenuation  measurements 
made  with  circular  piezoelectric  transducers  having  resonance  frequencies 
between  3  and  7  MHz. 
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A  comparison  Is  made  between  the  present  model  and  that  of 
Truell  and  Oates  [J.  Acoust.  Soc.  Am.  35,  1382  (1963)],  which  was 
Intended  as  an  Indicator  of  the  necessity  for  applying  a  correction 
for  nonparallelism.  The  present  model  corrects  the  first  four  echoes 

_3 

for  frequencies  f  <  6  MHz  and  angles  y  <  4  x  10  radians. 
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CHAPTER  I 

INTRODUCTION 

The  attenuation  of  ultrasonic  waves  is  a  physical  phenomenon 
which  has  great  significance  in  many  wave  propagation  studies.  Much 
useful  information  about  the  internal  structure  of  the  propagating 
medium  can  be  derived  from  a  correct  interpretation  of  attenuation. 

The  correct  interpretation  of  attenuation,  however,  is  a  long-standing 
problem.  The  usual  approach  to  solution  of  this  problem  is  to  assume 
that  the  physical  phenomena  which  produce  attenuation  make  independent 
contributions  to  the  measured  result.  One  then  characterizes  the 
measured  attenuation  coefficient  as  a  summation  of  coefficients  result¬ 
ing  from  (1)  intrinsic  attenuation  (that  resulting  from  viscous  and 
thermal  losses  in  the  propagating  medium  as  well  as  grain  boundary 
scattering);  (2)  attenuation  resulting  from  diffraction;  and  (3)  attenua 
tion  resulting  from  a  lack  of  parallelism  between  the  ultrasonic  source 
and  the  receiver.  This  thesis  Is  an  attempt  not  only  to  point  out  that 
(2)  and  (3)  are  not  truly  independent,  but  also  to  demonstrate  how  they 
can  be  interpreted  simultaneously. 

In  his  Theory  of  Sound  first  published  in  1878,  Lord  Rayleigh 
(1950)  gave  the  mathematical  form  of  the  diffraction  integral.  In 
1886  Lommel  published  a  solution  of  the  same  diffraction  integral  as 
applied  to  the  propagation  of  light  by  utilizing  the  Kirchhoff  approxi¬ 
mation  for  Fresnel  diffraction.  The  solution,  given  In  the  form  of  an 
infinite  sum  of  Bessel  functions,  is  awkward  to  use  (Gray  and 
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Matthews  (1931)).  Huntington  et  al.  (1948)  published  the  first 
diffraction  correction  for  propagation  of  ultrasonic  waves  by  making  a 
numerical  Integration  of  the  tabulated  data  of  Lommel.  Numerous 
authors  subsequently  have  provided  Improvements  of  the  numerical  Inte¬ 
gration  process,  the  latest  having  been  published  by  Benson  and 
Kiyohara  (1974).  Rogers  and  Van  Buren  (1974)  obtained  an  analytical 
solution  of  Lommel's  diffraction  integral.  Their  solution  is  applicable 
to  the  transmission  of  sound  between  two  discs  of  the  same  size  which 
are  parallel  to  each  other.  However,  the  analytical  solution  has 
limited  utility  in  the  present  situation  because  it  gives  the  result 
of  integrating  over  a  receiver  area  rather  than  the  value  of  pressure 
at  a  single  point  in  the  diffraction  field. 

All  the  above  models  assume  that  the  propagation  medium  is  a 
lossless  liquid.  Papadakis  (1975)  studied  the  propagation  of  ultra¬ 
sonic  waves  in  a  solid  and  states:  "The  solution  for  fluids  applies 
adequately  to  isotropic  solids  as  long  as  the  transducer  is  bonded 
adequately  to  one  of  a  pair  of  plane  parallel  faces  of  a  slab  consider¬ 
ably  larger  in  lateral  extent  than  the  transducer  diameter"  (p.  154). 

In  all  these  studies  the  receiver  and  the  source  conveniently  were 
assumed  to  be  parallel.  It  has  been  known  for  some  time  that  measure¬ 
ments  of  ultrasonic  attenuation  are  seriously  affected  by  nonparallelism 
of  the  source  and  the  receiver.  In  their  review  article  Breazeale  et  al. 
(1981)  state  "among  the  less  easily  evaluated  errors  are  those  arising 
from  .  .  .  phase  cancellation  resulting  from  lack  of  parallelism  of  the 
transducers,  from  material  Inhomogeneity,  and  from  diffraction"  (p.  68). 
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Truell  and  Oates  (1963)  studied  the  effect  of  nonparallelism  of  source 
and  receiver.  Their  study  showed  that  the  diffraction  pattern  observed 
in  a  sample  of  wedge  angle  y  was  modulated  after  n  reflections  by  a 
function  of  the  form  2Jj((2ka  ny)/2ka  ny).  They  then  gave  a  mathematical 
treatment  to  explain  this  modulation  phenomenon  on  the  basis  of  phase 
cancellation  of  infinite  plane  waves.  Truell  et  al.  (1969)  suggested 
that  the  form  of  the  function  could  be  used  as  a  means  of  predicting  the 
degree  of  parallelism  required  for  a  given  value  of  ka  to  avoid  the 
modulation  effect.  They  also  showed  that  the  model’s  ability  to  predict 
the  shape  of  the  modulation  was  much  better  at  high  frequency  (85  MHz) 
than  at  low  frequency  (35  MHz).  Calder  (1978)  introduced  a  modification 
to  the  Truell  and  Oates  model  by  correcting  for  the  inhomogeneity  of 
the  sample.  He  applied  his  correction  to  a  wedged  sample  and  to  a 
hemispherical  reflector.  In  his  case  the  Truell  and  Oates  model 
significantly  corrected  for  the  modulation  in  the  case  of  the  hemi¬ 
spherical  reflection;  however,  it  was  not  able  to  correct  for  the 
modulation  observed  with  the  wedged  sample. 

In  this  thesis  the  Truell  and  Oates  model  is  shown  to  be 
inappropriate  for  the  frequency  range  3  to  9  MHz  in  wedged  steel  samples. 
Next,  the  Truell  and  Oates  model  is  modified  by  the  removal  of  the 
plane  wave  approximation.  This  procedure  makes  it  necessary  to  define 
the  pressure  at  a  point  In  the  diffraction  field.  Then,  the  point 
values  of  the  pressure  are  numerically  Integrated  over  the  receiver 
area.  Such  an  integration  can  be  performed  only  after  a  definition 
of  the  relative  position  of  the  receiver  and  the  displaced  and  rotated 
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wave  field  has  been  made.  An  Inherent  characteristic  of  this  procedure 
Is  the  lack  of  separability  of  the  effects  of  diffraction  and  non¬ 
parallelism  of  the  sample  surfaces.  The  result  is  an  improved  model 
which  is  able  to  account  for  a  number  of  experimental  observations. 


CHAPTER  II 


THEORETICAL  ANALYSIS 
A.  WAVE  PROPAGATION  THEORY 


1.  Wave  Attenuation 

a.  Definition  of  a,  the  attenuation  parameter.  If  a  pressure 
disturbance  of  the  form  P(t)  *  pQe^  propagates  a  distance  Ax  from  a 
point  Xq  to  a  point  x,  then  the  pressure  measured  at  the  point  x  will 
be  of  the  form: 

P(x,t)  *  p0ej(&,t‘kAx)  .  (2-1) 

An  expression  for  an  attenuated  wave  is  obtained  by  assuming  that 
the  wave  vector  k  is  complex  and  given  by: 


k  = 


(2-2) 


where  a  is  the  attenuation  coefficient  which  is  assumed  to  be  constant. 
Equation  (2-1)  then  becomes 


„AY  j(wt-knAx) 
P(x,t)  «  [e  aAx]  PQe  C 


(2-3) 


The  bracketed  term  expresses  the  manner  In  which  the  attenuation 
coefficient  a  affects  the  amplitude  values.  If  we  now  take  the  modulus 
of  Eq.  (2-3),  we  obtain 
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pmxM 


P  e"aAx 
KOe 


(2-4) 


From  Eq.  (2-4),  then,  a  can  be  expressed  as 


«  =  -sr’n 


Wx> 


L  V 


or  equivalently  as 


-  20  Inn 

“  '  a5T  log10 


(2-5) 


(2-6) 


The  unit  of  a  in  this  case  is  db  per  unit  of  length  of  Ax, 

b.  The  attenuation  correction.  If  one  wants  to  correct  for 
attenuation,  he  must  first  examine  the  physical  causes  of  this 
phenomenon.  In  general,  the  causes  of  attenuation  are  divided  into  two 
categories:  the  intrinsic  attenuation  phenomena  and  the  geometric 
attenuation  phenomena  ag.  These  two  categories  are  assumed  to  be 
Independent  of  each  other  such  that  we  can  express  a  as: 


“  *  °1  +  °g  *  (2*7) 

The  intrinsic  attenuation  Is  a  measure  of  the  loss  of  energy 
of  the  pressure  wave.  In  a  solid  which  Is  reacting  as  a  continuum  most 
of  the  intrinsic  attenuation  measured  comes  from  the  phase  delay  between 
the  density  function  and  the  pressure  function.  Therefore,  In  a  given 
homogeneous  isotropic  medium  the  total  intrinsic  attenuation 
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experienced  by  a  wave  propagating  at  a  frequency  fQ  can  be  expressed, 
at  least  In  the  linear  approximation,  as  a  constant  a^.  Then  to  obtain 
a  value  of  Intrinsic  attenuation  one  can  first  correct  a  measured 
pressure  value  for  geometrical  attenuation  and  use  this  corrected 
pressure  In  Eq.  (2-6). 

The  geometrical  attenuation  phenomenon  is  the  one  that  will  be 
studied  in  greater  detail  as  It  Is  significantly  affected  by  many 
geometrical  factors  such  as  parallelism  of  the  faces  of  the  solid  and 
diffraction  which  have  been  adjusted  during  the  course  of  this 
investi gati on. 

2.  The  Diffraction  Theory 

a.  The  diffraction  model.  We  will  assume  that  a  source  of 
pressure  disturbance  behaves  as  a  circular  piston  in  an  Infinite  rigid 
baffle.  Then  the  magnitude  of  the  disturbance  at  any  point  on  the 
source  can  be  described  either  In  terms  of  maximum  pressure  PQ  or 
maximum  velocity  Uq  as: 


P(t)  *  P0eJu>t 

(2-8) 

U(t)  *  jUQejut 

(2-9) 

where  u  Is  the  angular  frequency  of  the  disturbance.  If  this 
disturbance  is  now  allowed  to  propagate  Into  a  lossless  homogeneous 
Isotropic  fluid.  Its  radiation  pattern  can  be  described  by  the  use  of 
classical  wave  theory.  Using  Huygens'  principle  we  then  approximate 
the  piston  by  subdividing  the  area  into  small  point  sources  and  use  the 
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additivity  theorem  to  obtain  the  resulting  effect  of  all  these  small 
sources  at  a  field  point. 

Klnsler  et  al.  (1982)  give  the  radiation  pattern  of  a  small  source 
(of  any  shaped  which  is  located  in  an  infinite  baffle  as: 


dP(R,t) 


5 


ds.0  ej(wt-kR)j 


(2-10) 


where 

R  Is  the  distance  from  the  center  of  the  source  to  the  field 
point; 

pg  is  the  medium  density; 

C  is  the  velocity  of  propagation  of  the  pulsation  through  the 
medium; 

Ug  is  the  maximum  velocity  of  the  surface; 

k  Is  the  wave  vector; 

dP  is  the  pressure  value;  and 

dS'  Is  the  pulsating  surface. 

The  pressure  at  a  field  point,  defined  by  the  spherical  coordinate 
(Rg,e)  shown  in  Figure  II— 1 „  Is  given  by  the  integration  of  Eq.  (2-10) 
over  S', the  source  surface.  One  then  obtains  the  diffraction  Integral: 


P(R0,e,t) 


J  p0  CU0  [ 

Js. - R - dS 


(2-11) 


^Appendix  A  shows  why  for  small  sources  the  shape  does  not  make 
any  difference  to  the  radiation  pattern. 
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Figure  II-l.  Geometry  Used  In  Deriving  the  Radiation 
Characteristics  of  a  Piston  Source. 
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This  Integral  also  is  referred  to  as  Rayleigh's  integral  and  it  does  not 
have  a  closed  form  Solution.  However,  a  number  of  numerical  integrations 
of  this  integral  have  been  made  and  the  result  of  one  of  them  is  shown 
in  Figure  II-2. 

b.  The  pressure  on  the  axis.  The  pressure  in  the  field  along 
the  Z  axis  is  one  of  the  few  quantities  which  can  be  analytically 
obtained  from  the  piston  model  without  further  approximation.  From 
Eq.  (2-11): 


P(z,0,t) 


j  p0  CU0k  ej(o>t) 

T 


i; 


0  (r2  +  z2) 


TJJT 


since: 


(2-12) 


-Jk(zV)'/2  .  .  0 

3r  JE 


(77^ e 


(2-13) 


Then  we  obtain 


P(z,0,t)  =  pq  CUQ  eJ[e~J  -  e 


j“>tr  -jkz  _  -jk(z2+a2)1/2] 


(2-14) 


The  magnitude  of  this  expression  is  the  pressure  on  the  axis  Paxj$: 


Paxis  *  2p0  cuols1n4  kz[{1  +  (a/z)2)V2  -  H>l  •  (2-15) 

When  the  value  of  z/a  Is  much  larger  than  ka,  then  Pax^s  simplifies  to 
an  asymptotic  expression  Psph  which  Is  the  spherically  diverging  wave 
expected  at  large  distances: 
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pO  2 

sph  '  2z  cu0  ka 

(2-16) 

sph(z>t)  '  57  CU0  ■ 

(2-17) 

A  comparison  of  the  values  of  P  obtained  for  Eqs.  (2-15)  and  (2-17), 
shown  in  Figure  I I -3,  shows  that  P$p^  does  not  converge  to  within  10% 
of  Pax^s  until  a  distance  of  twice  the  farfield  distance  zf  is  achieved. 
The  parameter  z^  is  defined  as  the  axial  distance  to  the  last  maximum 
of  the  axial  pressure  and  is  given  by  Krautkramer  (1977)  as: 

zf  =  (4a2  +  x2)/4x  (2-18) 

where  a  is  the  radius  of  the  source  element  and  x  is  the  wavelength  of 
the  wave. 


c.  The  pressure  in  the  farfield.  An  analytical  form  of  the 
diffraction  integral  can  be  obtained  by  noting  that  for  RQ  »  a  the 
value  of  R  is  adequately  approximated  by:  (1)  Rg  in  the  denominator 
and  (2)  by  Rg  -  r  sinecosi|>  In  the  more  sensitive  phase  term.  Then 
Eq.  (2-11)  can  be  integrated  to  give: 


P(Rg.e.t) 


*0  CU0 


p  j (wt-kRg) 
ka^  e  u 


2Jj(ka  sine) 
ka  sine 


Then  from  Eq.  (2-16) 


P(Rg,e,t)  *  psph^Ro^ 


2J.|(ka  sine) 
ka  sine 


J (»t-kR0) 


(2-19) 


(2-20) 
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where  the  spherical  dependence  *s  expressed  by  Eq.  (2-11) 

when  Z,  the  axial  distance,  is  equal  to  Rq.  The  remaining  term  includes 
all  the  angle  dependence  and  is  commonly  referred  to  as  the  directivity 
function.  A  plot  of  [2J^(ka  sine)/(ka  sine)]  is  shown  in  Figure  I 1-4.  It 
is  important  to  note  that  in  the  farfield  the  angle  e  gives  the  approxi¬ 
mate  direction  of  propagation  of  the  wave  on  the  wavefront,  i.e.,  along 
the  radius  RQ  and  that  the  function  H(0j)  becomes  much  more  sensitive 
to  variations  of  ei  as  the  value  of  the  parameter  ka  is  increased.  The 
polar  diagram  obtained  from  Krautkramer  (1977)  and  included  in 
Figure  I 1-5  clearly  shows  this  relationship. 

d.  Average  pressure  on  the  receiver.  In  a  pulse  echo  system 
one  can  obtain  the  average  pressure  on  the  receiver  by  integrating 
Eq.  (2-11)  over  the  receiver  area.  Rogers  and  Van  Buren  (1974)  have 
solved  analytically  this  problem  using  the  Fresnel1  diffraction 
approximation  and  obtained  a  diffraction  correction  parameter  D  given 
by: 


D=l.  e'j(27t/S)[J0(2tr/S)  +  j  J-| (2ir/S)] 


(2-21) 


2 

where  S  =  2itZ/ka  is 


the  normalized  propagation  distance. 


hn  the  Fresnel  approximation  one  replaces  R  in  Eq.  (2-11)  by  Z 
in  the  denominator  and  by  Z  +  (r*  +  r'z  -  2rr'cos  4»')/2Z  where  the 
prime  coordinates  are  the  source  coordinates.  (See  Gray  and  Mathews 
(1931).) 


80% 

(a)  For  ka  =  5u  (f  =  3.6  MHz,  a  *  1.27  cm;  in  steel) 


100% 
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Figure  II-5.  Directional  Characteristic  H(e)  (from  Krautkramer 
(1977)).  ' 
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B.  GEOMETRICAL  CORRECTION  MODELS 

1 .  The  Plane  Wave  Correction  Model 

As  will  be  shown  in  the  results,  significant  deviations  from  the 
diffraction  model  are  observed  in  the  echo  train  of  a  wedged  sample. 
Truell  and  Oates  (1963)  suggested  a  model  to  explain  the  deviation  from 
the  diffraction  pattern  based  on  the  geometry  of  wedged  samples.  This 
model  will  be  described  here. 

a.  The  geometry  of  propagation  in  a  wedged  sample.  Let  us 
assume  that  a  sinusoidal  disturbance  of  frequency  fQ  originating  at  a 
time  tg  and  having  radial  symmetry  about  an  origin  0,  begins  propagating 
along  the  normal  to  the  upper  surface  of  a  wedge  of  angle  y  as  shown  in 
Figure  I I -6.  At  any  time  t  >  tQ  the  center  of  symmetry  of  this  disturb¬ 
ance  is  located  along  the  axis  of  propagation  of  this  wave;  i.e.,  at  the 
center  of  the  wavefront.  Now,  consider  the  displacement  of  the  center 
of  the  wavefront  of  the  disturbance  as  the  wave  is  multiply  reflected 
within  the  sample.  This  displacement  occurs  along  a  path  described  by 
a  ray  folded  by  reflection  Detween  the  faces  of  the  wedge.  At  every 
reflection  from  the  lower  face  the  relative  angle  vn  between  the  upper 
surface  and  the  normal  to  the  propagation  direction  of  the  wave  is 
Increased  by  2y.  After  n  reflections  from  the  lower  surface,  the 
relative  angle  Yn  is  given  by: 

Yn  *  2ny.  (2-22) 


Furthermore,  it  is  evident  that  as  the  wave  propagates  its  symmetry 
axis  Z  will  wander  away  from  the  receiver  origin  0. 
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b.  The  theory  of  Truell  and  Oates.  Truell  and  Oates  (1963) 
gave  a  theoretical  model  to  explain  the  rapid  fluctuation  of  amplitude 
observed  in  the  echo  train  of  a  circular  piezoelectric  receiver  when  a 
wedged  sample  is  used.  They  assumed  that  the  wavefront  was  effectively 
plane  and  that  the  fluctuations  in  pulse  height  could  be  accounted  for 
by  the  variations  in  the  receiver  signal  that  resulted  from  phase 
cancellation  of  the  wavefronts  across  the  receiving  transducer. 

Assuming  the  effective  translation  of  the  receiver  is  negligible,  then 
the  pressure  function  at  a  point  (Rg,0,t)  after  n  reflections  from  the 
lower  surface  of  the  wedge  can  be  expressed  as: 


j{wt-k[2nL-xsinY  ] 

P(R0.e,t)n  *  e 


(2-23) 


This  expression  must  be  integrated  over  the  circular  receiver  to  obtain 
the  average  pressure  (P)n: 


(P) 


n 


?  1/9  -jk(2nl-xsiny  ) 
xZ)1/Z  e  n  dx  . 


(2-24) 


The  result  of  this  integration  gives: 


(P)n  =  sin(oot  -  2knl)  H(yn) 

where 


(2-25) 


H(yn) 


aj^ka  s1nyn] 


ca  sinyr 


(2-26) 


The  form  of  this  function  is  given  by  the  magnitude  of  the  directivity 
function  defined  in  Eq.  (2-20). 


c.  The  average  pressure  on  the  receiver.  In  the  Truell  and 
Oates  model,  the  function  obtained  in  Eq.  (2-25)  is  assumed  to  modulate 
the  exponential  diffraction  pattern  such  that  both  effects  can  be 
corrected  for  separately.  The  expression  of  the  average  amplitude  on 
the  receiver  is  then  given  by  the  multiplication  of  the  diffraction 
correction  parameter  D  given  by  Eq.  (2-22)  and  H(yn).  Therefore,  the 
magnitude  of  the  average  pressure  P  is  given  by: 


|P|  -  DH(yn)  . 


(2-27) 


The  range  of  |P|  may  go  from  0  to  1  and  conversely  the  geometrical 
attenuation  according  to  the  plane  wave  model  goes  from  -®  to  0. 


2.  The  Diffracted  Wave  Model 

As  will  be  seen  In  the  Results  section,  the  plane  wave  model 
which  assumes  Independence  of  the  diffraction  and  wedge  geometry 
correction  does  not  satisfactorily  describe  the  echo  train  observed. 

A  new  approach  combining  both  the  diffraction  theory  and  the  description 
of  the  phase  cancellation  will  be  developed  for  the  wedge  sample 
geometry. 
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a.  Model  of  the  physical  situation.  If  one  examines  the 
diffraction  model  given  In  Eq.  (2-11)  and  depicted  by  Figure  II-2 
(p.  11),  one  realizes  that  assuming  a  plane  wave  with  a  constant 
pressure  distribution  Is  a  very  crude  approximation  of  the  physical 
situation.  Furthermore,  this  description  Is  made  In  relation  to  a 
fixed  propagation  axis.  However,  as  Is  obvious  from  Figure  II-6,  the 
propagation  axis  Is  rotated  and  translated  in  relation  to  the  source/ 
receiver  element  as  the  wave  Is  multiply  reflected  in  the  sample.  An 
example  of  an  analogous  situation  for  a  unidirectional  propagation  path 
Is  shown  schematically  In  Figure  II-7.  Using  this  physical  model  we 
will  define  the  average  pressure  resulting  on  the  receiver  when  the 
propagation  distance  Is  larger  than  the  nearfield  distance.  First,  a 
modification  of  the  farfield  approximation  will  be  done  to  improve  its 
convergence.  Second,  a  coordinate  transformation  system  will  be 
defined  such  that  the  unidirectional  diffraction  theory  can  be  applied 
to  the  propagation  In  a  wedge  sample.  Third,  an  average  pressure  will 
be  obtained  by  numerical  Integration  of  the  pressure  function  over  the 
receiver  area. 

b.  Modification  of  the  solution  of  the  diffraction  Integral  in 
the  farfield.  A  study  of  the  convergence  of  the  farfield  approximation, 
given  in  Eq.  (2-20),  to  a  numerical  approximation  of  the  diffraction 
Integral  was  made  by  Rose  (1975).  He  concluded  that  the  directivity 
function  (2J^(u)/u)  seemed  to  converge  much  more  rapidly  to  the  shape 
of  the  numerical  values  than  the  value  of  P  ^  to  the  value  of  pressure 
on  axis  Pax^s»  In  order  to  obtain  a  more  rapid  convergence  of  the 
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solution,  the  value  of  P$ph  In  Eq.  (2-20)  is  replaced  by  the  value  of 
Paxis  *ron  {2-15).  The  new  equation  for  the  pressure  Is  then: 

j(ut-kRn) 

P(R0,e,t)  «  Paxis(R0)[2J1(kas1ne)/(kasine}]e  u  .  (2-28) 

A  numerical  integration  of  this  integral  was  done  over  a  receiver 
of  1.27  cm  radius  using  the  IMSL  function  DBLINT  which  is  a  cautious 
Romberg  quadratic  numerical  integration  process.  This  result  is  com¬ 
pared  with  the  closed  form  analytical  expression  of  the  diffraction 
Integral  given  by  Eq.  (2-22).  Figure  II-8  shows  the  value  of  Eq.  (2-28) 
as  a  means  of  improving  the  convergence  of  the  expression  of  the 
pressure  in  the  farfield. 

The  normalized  plot  of  average  pressure  shown  in  Figure  I I -8  is 
then  used  to  extend  the  validity  of  the  numerical  Integration  of 
Eq.  (2-28)  to  points  having  a  propagation  distance  smaller  than  three 
times  Zf.  A  correction  Is  then  applied  to  all  the  average  pressure  P 
subsequently  obtained  by  Integrating  this  equation  over  the  receiver 
surface.  The  correction  Is  applied  by  multiplying  the  average  pressure 
by  the  ratio  of  the  average  pressure  obtained  from  Rogers  and  Van  Buren's 
solution  to  the  average  pressure  obtained  by  integrating  Eq.  (2-28)  for 
a  piezoelectric  element  of  1.27  cm  radius  at  4  MHz.  This  radio  was 
verified  to  be  the  same  to  within  less  than  K  whether  the  radius  is 
0.635  cm  or  1.27  cm  or  the  frequency  2,  4,  6,  or  8  MHz. 

c.  The  coordinate  transformations. 

The  new  coordinate  system.  To  define  a  new  coordinate  system 
it  will  be  necessary  to  define  an  imaginary  source  of  origin  0'  located 


Figure  I 1-8.  Inverse  of  the  Average  Pressure  on  the  Receiver 
Function  of  normalized  Propagation  Distance. 
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along  the  axis  of  symmetry  Z'  of  the  wave  such  that  an  equivalent 

unidirectional  path  can  be  associated  with  the  wave.  In  this  system 

the  axis  Z'  Is  defined  to  be  coincident  with  the  axis  of  symmetry  and 

the  direction  of  propagation  of  the  wave  after  n  reflections  from  the 

lower  face.  From  symmetry  consideration  It  will  then  be  mathematically 

simpler  to  use  a  cylindrical  coordinate  system  described  by  the  three 

orthogonal  axes  p',  and  Z'. 

Position  of  the  Imaginary  source  origin  0*.  The  radial 

coordinate  [p'(0,0)]n  of  the  center  0  of  the  receiving  element  is 

obtained  by  considering  the  right  triangle,  formed  by  the  intercept  of 

the  Z  axis,  the  Z'  axis  and  the  distance  [p'(0,0)]n.  In  Figure  I 1-9. 

The  angle  between  the  normal  to  the  surface  and  the  axis  Z'  is  defined 

by  Eq.  (2-22)  as  Yn«  Its  opposite  angle,  defined  by  the  axes  Z  and  Z', 

must  also  be  since  Z  is  normal  to  the  surface.  Then  we  have 
n 

[p-(0,0)Jn  ■  n  L  s1nyn  ,  (2-29) 

where  L  is  the  sample  length  measured  at  the  center  of  the  receiver. 

The  angle,  defined  by  the  distance  [p'(0,0)]n,  and  the  surface  must  be 
equal  to  yn  since  Its  two  sides  are  respectively  perpendicular  to  the 
normal  to  the  surface  and  to  the  Z'  axis.  Then  its  bisectrix  will  cut 
the  Z'  axis  at  a  point  [Z'arJn.  Consider  now  the  triangle  formed  by 

A 

the  bisectrix,  the  Z  and  the  Z'  axis.  In  this  triangle  the  angle  A 
opposite  to  the  Z'  axis  side  Is  given  by 
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Figure  I 1-9.  Geometrical  Relationship  Existing  for  n  *  1. 
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Since  the  angle  B  opposed  to  the  bisectrix  has  been  defined  previously 
as  Yn,  then  the  angle  C  opposite  to  the  Z  axis  is  given  by: 

C  •  »  -  B  -  A  (2-31) 

or 

c  -  IT  -  Yn  -  (*/2  -  Yn/2)  =  n/2  -  yn/2  .  (2-32) 

Since  C  equals  A,  then  their  sides  are  equal  and  since  the  Z  axis  side 
is  equal  to  n  x  L,  then  the  axis  side  is  also  equal  to  n  x  L.  This 
means  that  the  propagation  path  of  the  wave  from  0  to  [Z'ar]n  1S  given 
by: 


[Zpar^n  =nxL+nxL=2nL  •  (2-33) 

The  magnitude  of  [Z'ar]n  is  the  propagation  distance  of  a  wave  between 
parallel  faces  after  n  reflections  from  the  lower  surface. 

Now  let  us  consider  the  right  triangle  formed  by  the  distance 
[p'(0,0)]n,  the  bisectrix  and  the  Z^  axis.  The  Z'  axis  side  which  we 
will  call  hi',  is  given  by: 

AZ-  =  [p-(0,0)]n  tan(yn/2)  .  (2-34) 

The  propagation  path  from  0  to  [Z(0,0)]n  is  given  by: 

[Z'(0,0)]n  =  [Z'ar]n  -  hi'  .  (2-35) 

Then  by  combining  (2-33),  (2-34),  and  (2-35),  we  obtain  the  [Z'(0,0)]n 
coordinate  as: 
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[Z'(0,0)]n  =  n  x  L  [2  -  sin(yn)  tan(yn/2)]  .  (2-36) 

The  relative  positions  of  the  imaginary  source  origin  0'  and  the 
receiver  origin  0  are  now  defined  by  Eqs.  (2-29)  and  (2-36).  This 
geometrical  relationship  which  has  been  shown  to  be  exact  for  the  first 
reflection  is  assumed  to  be  exact  for  n  reflections.  It  has  been 
verified  graphically  to  be  exact  for  at  least  four  reflections  for 
different  angles  y  ranging  from  1°  to  5°.  An  example  of  such  graphical 
construction  is  shown  in  Figure  11-10. 

Cylindrical  coordinate  of  a  point  on  the  receiver  surface. 
Consider  the  projection  of  the  upper  surface  in  the  p'  Z'  plane  cut 
along  the  X  axis  as  shown  in  Figure  11-11.  The  X  axis,  the  parallel  to 
the  p'  axis  passing  by  the  poing  (0,0)  and  the  parallel  to  the  Z'  axis 
passing  by  the  point  (x,0)  form  a  right  triangle  of  respective  side  x, 
xcosyp  and  xsinyn.  Then  the  intercept  on  the  Z'  axis  of  the  point 
(x,0)  is  given  by 

[Z'(x,0)]n  =  [Z'(0,0)]n  -  xsinyn  .  (2-37) 

The  points  (x,y)  and  (x,0)  must  have  the  same  Z'  value  since  y 
is  antiparallel  to  Therefore,  Eq.  (2-37)  can  be  rewritten  as: 

[Z'(x,y)]n  =  [Z'(0,0)]n  -  xsinyn  .  (2-38) 

The  intercept  of  the  point  (x,0)  on  the  p'  axis  is  similarly  given  by: 

[p'(x,0)]  =  [p'(0,0)]n  +  xcosyn  . 


■  --  •  ♦  - «A3T.  ■  '  lm  »)'  'll*  .  ■  ■  Im  ,  ..  ,  in'  m 
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Figure  11-10.  Graphical  Construction  Showing  the  Geometrical 
Relations  after  Three  Reflections. 


X 


V" 


Figure  11-11.  Definition  of  the  Receiver  Coordinate  (x,y)  i 
the  o'Z'  Plane. 
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Let  us  now  consider  the  p'  plane  cut  along  the  parallel  to  the  y  axis 

formed  by  the  point  (x,0)  and  (x,y)  as  shown  in  Figure  11-12.  In  the 
right  triangle  formed  by  the  three  points  (x,0),  (x,y),  and  [0,Z'(x,y)]n, 
the  distance  between  the  last  two  points  is  [p'(x,y)]n  and  is  given  by: 


[p'(x.y)]n  =  CCp'(x,0)]2  +  y2]1/2 


We  can  now  convert  x  and  y  in  terms  of  r  and  41  by  noting 


x  =  r  cos  i|j  , 


(2-39) 


y  =  r  sin  ip  . 


(2-40) 


(2-41) 


This  means  that  we  can  define  the  diffraction  field  parameters  Rq,0  in 
terms  of  r  and  ij>  by  using  the  following  identities: 


(R0)  -  {[p'(r,*)l2  +  [Z'(r,*)]2}1/2  , 

.1  [p'(r,i}*)3 
/  ^  \  _  a. _ "  1  n 


(e)n  *  tan 


(2-42) 


(2.43) 


d.  The  average  pressure  on  the  receiver.  From  the  modified 
expression  of  the  pressure  P  at  a  point  in  the  field  expressed  by 
Eq.  (2-28)  and  the  relationship  between  Rq,  e  and  (p',Z')n  defined  by 
Eqs.  (2-42)  and  (2-43),  we  can  express  the  pressure  at  a  point  (p',Z')_ 
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P(o',z',t)n  =  Paxis(R0(p>,2')n)[2J1tkas1ne(p',z')n)/icasine(p',z')n] 
5  (wt-kR#j (p',z')  ) 

e  0  n  .  (2-44) 


By  using  Eqs.  (2-38)-(2-42)  we  obtain  p'  and  z'  as  a  function  of  r  and 
iji,  the  receiver  coordinates.  The  average  pressure  Pn(t)  is  obtained  by 
integrating  Eq.  (2-44)  over  the  circular  receiver  area  of  radius  a. 

This  gives: 


pn(t>  --tT  r 

n  ita  Ji>= 0  Jr=0 


j[wt-kR0(r,K-)] 


r  dr  dif/  . 


(2-45) 


Since  in  the  attenuation  study  we  measure  the  maximum  pressure,  we  need 
the  magnitude  of  Eq.  (2-45).  Using  the  identity 

eje  *  cose  +  j  sine  ,  (2-46) 

Eq.  (2-46)  becomes 

Pn(t)  *  "T  |0  lQ  lpCr,ip] lnfcos[u>t  -  kRQ(r,i|i)] 

+  j  sin[wt  -  kRg(r,if/)]  )r  dr  dxp  .  (2-47) 

Since  we  need  only  the  real  part,  then 

P„(t)  *  "V  f  [  I P[r, *]  (  cos [wt  -  kR0(r,i//)]r  dr  dif/  .  (2-48) 

n  na^  Jo  Jo  n  0 
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Since  cos(A-B)  Is  equal  to  cosA  cosB  +  sinA  sinB,  then 


~  i  fa  fa 

P_(t)  ■  — jr  (coswt  !P[r,^]LcoskR0(r,4i)r  dr 

n  va*  h  >0  n  u 


dip 


fa  f3 

sinwt  j  j  |P[r,i(i]jnsinkR0(r,if))r  dr  d<f>}  •  (2-49) 


These  two  integrals  can  now  be  evaluated  numerically  such  that 
Eq.  (2-49)  becomes 


Pp(t)  =  (iTaZ)"1{cosa)t  +  sinwt  fg)  •  (2-50) 

The  magnitude  of  Pn(t)  is  given  by 

|P|„  -  [.a2]"'  x  [f,2  a  f22]1/2  .  (2-51) 

According  to  the  diffracted  wave  model,  the  parameter  | P j n  then 
represents  the  maximum  average  pressure  impinging  on  a  circular  trans¬ 
mitter/receiver  when  a  pressure  wave  subjected  only  to  geometrical 
attenuation  and  of  a  magnitude  of  1  is  reflected  n  times  from  the  lower 
face  of  a  wedged  sample.  The  range  of  | P i n  goes  from  1  to  0  and  the 
corresponding  range  of  (the  attenuation  due  to  geometrical  attenua¬ 
tion)  goes  from  0  to  —>. 

3.  The  Correction  for  Geometrical  Attenuation 

An  amplitude  correction  factor  (CG)n  is  defined  as  the  ratio  of 
the  magnitude  of  the  average  pressure  felt  on  a  receiver  of  radius  a 
when  (1)  an  infinite  plane  wave  is  impinging  normally  on  its  surface 
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(P  *  1)  to  (2)  the  situation  described  by  either  Eq.  (2-51)  or 
(2-27). 


(CG)  *  -J-  .  (2-52) 

|P|„ 

Then  after  n  reflections  from  the  lower  face  of  a  wedge  the  magnitude 
of  a  measured  average  pressure  |P|n  can  be  corrected  for  geometrical 
effects  to  obtain  a  corrected  average  pressure  |Pcln- 


|PCI„  *  |P|„  X  I CG i 


(2-53) 


Equivalently,  |Pcln  is  the  magnitude  of  an  infinite  plane  wave 
propagating  normally  to  the  receiver.  The  pressure  P(t)  in  such  a  wave 
is  given  by: 


p(ti  -  ipci„ 


j(»t-kR0(0,0)n) 

e 


(2-54) 


where  Ro(°»0)n  is  the  distance  from  the  center  of  the  imaginary  source 
0'  to  the  center  of  the  receiver  (0,0). 

The  geometrical  attenuation  correction  process  described  In  this 
section  is  coded  in  FORTRAN  10  language  and  Included  in  Appendix  B. 

The  program  calculates  the  attenuation  between  the  first  and  the  nth 
corrected  echo  of  a  given  pulse  echo  setup  specified  by:  (1)  the 
velocity  V  of  the  ultrasonic  wave  in  the  sample,  (2)  the  source  radius 
a,  (3)  the  frequency  of  the  piezoelectric  crystal  f,  (4)  the  non¬ 
parallelism  angle  y»  (5)  the  sample  length  L,  and  (6)  the  relative 
amplitude  of  the  voltage  from  the  piezoelectric  crystal  after  n 
reflections  A(n). 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A.  PURPOSE  AND  BACKGROUND 

The  purpose  of  this  Investigation  is  to  gain  insight  into  the 
propagation  of  ultrasonic  waves  in  a  wedged  sample.  Preparatory  efforts 
were  concentrated  in  three  areas:  sample  definition  and  characteriza¬ 
tion;  design  of  electronic  setup  for  transmission  and  reception  of 
electrical  signal;  design  of  a  transducer  for  coupling  of  ultrasonic 
waves  and  electrical  signal. 

B.  EXPERIMENTAL  APPARATUS 
1.  Electronic  Components 

a.  System  description.  The  block  diagram  of  the  electronic 
system  is  shown  in  Figure  III.l.  A  variable  frequency  oscillator  (VFO) 
is  used  to  generate  an  RF  signal  of  a  given  frequency.  This  signal  is 
then  fed  into  a  frequency  counter  and  a  gated  pulsed  amplifier.  The 
pulsed  amplifier  sends  out  amplified  pulses  of  this  RF  signal,  at  a 
specified  pulse  repetition  frequency  (PRF),  to  the  diode  expander.  The 
diode  expander  eliminates  undesirable  low  amplitude  secondary  pulses  by 
the  use  of  Zener  diodes.  The  signal  is  then  fed  into  a  decoupler 
limiter.  This  decoupler  sends  the  signal  to  the  receiver  and  to  the 
transducer  and  limits  the  amplitude  of  the  signal  fed  to  the  receiver 
to  avoid  saturation  of  the  receiver  amplifier.  The  transducer 
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Figure  I I 1.1.  Electronic  Circuits  Used  for  Acoustical  Measurements. 
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transforms  these  single  frequency  pulses  into  longitudinal  pressure 
waves  of  ultrasonic  frequency  corresponding  to  the  frequency  sent  by 
the  VFO.  As  the  wave  propagates,  an  echo  is  eventually  reflected  back 
onto  the  transducer  where  the  piezoelectric  element  transforms  this 
pressure  wave  back  into  an  electrical  output  having  the  frequency  of  the 
original  signal,  but  of  much  lower  amplitude.  This  echo  is  then  retrans¬ 
mitted  by  the  decoupler  limiter  to  the  receiver.  The  receiver  contains 
a  wide  band  (1  to  90  MHz)  amplifier  which  can  amplify  received  RF  signals 
by  up  to  20  decibels.  This  signal  is  then  rectified  and  can  be  amplified 
by  up  to  70  decibels  by  the  second  stage  of  this  receiver.  The  rectified 
and  amplified  transducer  response  is  then  fed  into  the  boxcar  integrator 
and  displayed  on  the  oscilloscope.  The  position  of  that  gate  relative 
to  the  signal  indicates  the  part  of  the  signal  that  is  averaged  out  and 
measured  by  the  boxcar  integrator.  The  synchronization  of  the  pulser- 
receiver-boxcar  integrator  system  is  ensured  by  triggering  all  these 
elements  from  the  pulser  internal  trigger.  The  oscilloscope  is  then 
triggered  from  the  boxcar  delayed  trigger  output. 

b.  System  discussion.  Pulse  operation  was  selected  over 
continuous  wave  for  its  superior  immunity  to  cross-talk  interference. 
Pulse-echo  was  selected  over  through  transmission  to  minimize  system 
alignment  error.  The  need  for  accurate  frequency  measurements  and 
short  pulse  length  suggested  the  use  of  a  VFO  coupled  to  a  gated 
amplifier  rather  than  an  integrated  oscillator-pulser-receiver  unit. 

As  mentioned  by  Truell  et  al.  (1969),  this  type  of  setup  allows  for 
excellent  sensitivity  to  change  in  amplitudes  for  all  values  of 
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attenuation.  A  limiter  Is  built  Into  the  receiver  Input;  however,  to 
avoid  partial  saturation  of  the  receiver  amplifier  it  is  still  necessary 
to  use  an  external  decoupler  limiter.  A  diode  expander  is  used  to 
eliminate  pulses  of  low  level  energy  which  appear  Immediately  after  the 
collapse  of  the  main  RF  pulses1  and  last  long  enough  to  affect  the 
first  two  returning  echoes.  Terminators  of  93n  impedance  are  used  to 
provide  impedance  matching  between  the  pulser  and  the  transducer. 

2.  Transducer 

a.  Structure.  A  quartz,  x-cut,  crystal  is  used  as  a 
piezoelectric  element  for  the  generation  and  reception  of  the  ultrasonic 
wave.  As  discussed  in  Krautkramer  (1969),  one  of  the  chief  advantages 
of  quartz  crystal  over  ceramic  crystals  such  as  barium  titanate  or 
lithium  niobate  is  its  mechanical  strength  which  was  necessary  to  allow 
for  pressure  application  of  the  crystal  on  the  sample.  This  pressure 
is  applied  by  the  use  of  a  vise  pressing  on  both  the  sample  lower  face 
and  the  transducer  brass  housing  at  the  same  time,  as  can  be  seen  in 
Figure  I I I. 2.  A  photograph  of  the  transducer  parts  is  also  shown  in 
Figure  III. 3.  The  transducer  parts  were  originally  designed  to  be  used 
with  different  types  of  Immersion  transducers  [see  Scott  (1975)].  The 
main  brass  housing  is  maintained  by  the  vise  on  a  lucite  crystal  housing. 
This  crystal  housing  is  of  the  form  of  a  ring  with  the  upper  side  flat 
against  the  brass  housing  and  a  depression  on  the  lower  side.  This 

Hhls  effect  Is  generally  caused  by  accumulation  of  low  level 
piezoelectric  energy  In  the  ceramic  capacitor  of  the  oscillator  circuit, 
which  is  coupled  to  the  load  after  the  collapse  of  the  RF  pulse. 
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Figure  I I 1-2.  Ultrasonic  Measurement  Setup  Showing  the  Steel 
Plate,  the  Transducer,  and  the  Vise. 


41 


Figure  I I 1-3-  Transducer  Parts.  Top  row,  left  to  right;  lucite 
crystal  housing  ring,  lucite  cylinder  housing,  contact  spring,  contact 
cylinder.  Bottom  row:  main  brass  housing. 
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depression  has  an  inside  radius  1  mm  smaller  than  the  crystal  radius  and 
an  outside  radius  of  about  the  size  of  the  crystal  so  that  the  crystal 
can  be  inserted  in  this  depression.  Since  the  depth  of  this  depression 
is  smaller  than  the  thickness  of  the  crystal,  the  pressure  applied  on 
the  housing  is  uniformly  distributed  on  a  1  mm  ring  at  the  rim^  of  the 
crystal.  It  was  necessary  to  use  pressure  and  wring  the  crystal  on  the 
surface  for  two  main  reasons.  First,  it  was  a  simple  yet  efficient  way 
to  ensure  parallelism  between  crystal  and  sample  faces.  As  pointed  out 
by  Truell  et  al.  (1969).,  "ordinary  bonding  techniques  do  not  necessarily 
provide  sufficiently  good  control  over  parallelness"  (p.  120).  Further¬ 
more,  the  effect  of  nonparallelism  of  sample  faces  is  almost  impossible 
to  differentiate  from  the  effect  due  to  lack  of  parallelism  between  the 
crystal  and  the  sample  faces.  Secondly,  to  wring  the  crystal  on  the 
surface  is  a  very  efficient  way  to  damp  the  crystal  and  avoid  the  "ring¬ 
ing"  effect  described  by  Krautkramer  (1977).  In  this  case  the  cohesive 
force  between  the  sample  and  the  crystal  allows  the  sample  itself  to  act 
as  damping  material.  This  negates  the  need  for  rubber  backing,  the 
usual  damping  material,  which  causes  a  significant  decrease  of  crystal 
sensitivity. 

b.  Electrical  circuit.  The  electrical  signal  from  the  decoupler 
enters  at  the  BNC  connector  of  the  main  housing.  It  is  then  fed  through 
the  central  spring  loaded  retractable  connector.  This  connector  is  in 

^y  applying  the  pressure  at  the  rim  of  the  crystal  rather  than 
at  the  center,  uneven  loading  of  the  crystal  was  avoided  as  well  as 
unnecessary  loss  of  sensitivity  to  returning  echoes. 
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contact  with  a  brass  cylinder.  This  cylinder  can  move  up  and  down 
within  the  lucite  cylinder  housing  which  fits  into  the  center  of  the 
crystal  housing  ring.  The  cylinder  is  spring  loaded  by  a  central 
connector,  located  inside  the  brass  housing,  so  that  it  stays  in  con¬ 
tact  with  the  electrode  attached  to  the  surface  of  the  crystal. 

Two  types  of  electrodes  are  used:  gold  and  silver J  The  electrical 
circuit  is  then  completed  by  connecting  the  lower  face  of  the  crystal 
to  the  brass  housing.  This  is  achieved  by  connecting  the  sample  to  the 
brass  housing  using  a  brass  connector. 

3.  Sample 

The  sample  used  is  a  plate  of  high  carbon  steel.  This  material 
was  selected  for  its  known  characteristics  of:  isotropy,  low  attenuation 
coefficient  from  2  to  5  MHz,  electrical  conductivity  and  machinability. 
The  plate  geometry  is  necessary  to  eliminate  side  effects  which  can 
easily  be  confused  with  nonparallelism  effects  as  noted  by  Truell  et  al . 
(1969),  and  to  facilitate  machinability  of  specified  angles. 

The  plate  was  first  ground  to  ensure  parallelism  of  upper  and 

_2 

lower  faces  and  to  eliminate  surface  defects  larger  than  1  x  10  cm 
(this  represents  about  one-tenth  of  the  wavelength  at  6  MHz).  The  upper 
plate  was  then  ground  in  facets  of  approximately  3  cm  in  width.  A 
drawing  showing  the  approximate  geometry  and  the  coordinate  system  used 
is  shown  in  Figure  III. 4.  The  z  direction  corresponds  approximately  to 

^  gold  electrode  was  plated  on  the  entire  (1.27  cm-radius)  crys¬ 
tal  surfaces.  For  the  7  MHz  crystal  only,  a  1.14-cm  partial  aluminum 
electrode  was  attached  by  a  grease  couplant  which  was  squeezed  as  thin 
as  possible.  Such  partially-plated  transducers  have  been  studied  by 
Papadakis  (1975). 
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the  propagation  direction  of  the  echoes,  the  y  direction  corresponds  to 
the  width  of  the  sample  and  the  -x  direction  corresponds  approximately 
to  the  direction  of  lateral  displacement  of  the  propagating  wave.  As 
we  go  along  the  x  direction,  we  note  an  increase  of  the  relative  angle 
between  the  upper  facets  and  the  lower  face.  This  relative  wedge  angle 
(y)  increases  from  0  to  45  minutes  by  increments  ranging  from  20  seconds 
to  about  15  minutes.  The  width  Wn  of  the  facets  ranges  from  3  to  5  cm. 

C.  EXPERIMENTAL  PROCEDURE 
1 .  Sample  Characterization 

a.  Angles  measurements.  The  angles  between  the  lower  face  and 
each  facet  is  evaluated  by  using  a  collimator  and  two  reflectors  as 
shown  in  Figure  I I I. 5.  The  collimator  and  the  first  reflector  are 
positioned  such  that  the  cross  hairs'  image  coming  from  the  lower  face 
reflector  is  coincident  with  the  micrometer  cross  hairs.  The  micrometer 
cross  hairs'  position  is  noted.  The  second  reflector  is  now  positioned 
at  the  edge  of  the  facet  of  the  upper  face  to  be  measured  and  the 
micrometer  cross  hair  positioned  to  be  coincident  with  the  image  from 
this  upper  face  reflector.  Since  both  reflectors  have  parallel  and  flat 
faces  and  have  one  or  part  of  one  of  their  faces  wrung  onto  the  sample, 
the  wedge  angle  (y)  between  both  faces  is  given  by: 
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Lower  Face  Reflector 


Figure  III. 5.  Schematic  of  Optical  Measurement  Setup. 
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where 

y  *  wedge  angle; 

f  =  focal  length^  of  the  collimator  ■  33.23  mm; 
tx  =  distance  between  the  two  image  positions  along  the  x  axis. 
The  distance  Ln  is  the  sample  length  measured  between  the  lower 
end  and  the  upper  facet  n.  Four  measurements  were  made  at  each  facet: 
Yn  at  both  edges  of  the  sample  (y  *  0  cm  and  y  =  10.5  cm)  and  Ln  at 
about  3  cm  from  either  edge  of  the  sample  (y  =  3  cm  and  y  =  7  cm).  The 
position  of  y  for  measurements  of  Ln  correspond  approximately  to  the 
position  of  the  center  of  the  transducer.  The  result  of  these  measure¬ 
ments  Is  shown  in  Table  III.l. 

b.  Velocity  measurement.  The  measurement  of  velocity  is  made 
by  using  a  method  analogous  to  the  long  pulse  technique  described  by 

2 

McSkimin  (1950).  First,  the  facet  used  for  measurements  is  thoroughly 
cleaned  with  acetone.  The  crystal  is  then  positioned  at  the  center  of 

3 

this  facet  on  a  drop  of  couplant  squeezed  between  the  crystal  and  the 
sample  to  ensure  a  good  mechanical  coupling.  The  transducer  housing 


Hhe  value  of  f,  the  focal  length,  was  measured  using  the 
Newtonian  method  described  by  Palmer  (1969). 

2 

The  cleaning  Is  extremely  Important,  as  it  prevents  lack  of 
parallelism  or  bad  coupling  due  to  the  presence  of  dust  particles. 

3 

After  many  trials,  nonaq  stopcock  grease  was  selected  as 
couplant.  Its  chief  advantage  Is  that  when  pressure  was  applied,  an 
even  and  very  thin  layer  of  couplant  was  formed.  This  permitted  to 
wring  the  crystal  very  easily  on  the  sample  and  the  ultrasonic  trans¬ 
mission  properties  of  this  coupling  were  excellent.  Another  advantage 
Is  that  the  crystal  could  afterwards  be  easily  removed  by  sliding  it  to 
the  edge  of  the  sample. 
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TABLE  III. 1 
SAMPLE  GEOMETRY  DATA9 


n 

Yn  (mil  11  radians)*5 

Ln  (cm. )c 

y  =  10.0  cm 

y  =  3.0  cm 

y  =  7.0  cm 

1 

0.14 

0.14 

2.4845 

2.4845 

2 

0.11 

0.19 

2.4855 

2.4855 

3 

0.00 

0.00 

2.4856 

2.4865 

4 

0.27 

0.16 

2.4860 

2.4865 

5 

0.28 

0.00 

2.4860 

2.4865 

6 

1.60 

1.62 

2.4845 

2.4850 

7 

3.82 

3.87 

2.4773 

2.4775 

8 

8.75 

8.07 

2.4563 

2.4570 

9 

11.64 

11.75 

2.4243 

2.4280 

aThe  variations  of  yn  and  Ln  observed  along  the  y  direction  were 
caused  by  machining  problems. 

h  .0 

Accuracy  of  the  measurements,  ±  3  x  10  mill i radians. 

c  -5 

Accuracy  of  the  measurements,  ±  5  x  10  cm. 
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is  then  mounted  on  the  crystal  and  the  pressure  slowly  increased  by 
tightening  the  vise.  During  this  operation  the  RF  pulse  rather  than  the 
boxcar  integrator  gate  position  is  monitored  on  the  oscilloscope.  (This 
connection  Is  shown  as  a  dotted  line  on  the  block  diagram  of  Figure  III-l, 
p.  37.)  When  the  RF1  and  rectified  pattern,  as  seen  on  the  oscilloscope, 
show  no  signs  of  distortion  and  when  increasing  the  pressure  does  not 
improve  the  pattern,  the  boxcar  integrator  gate  position  indicator  out¬ 
put  is  reconnected  to  the  oscilloscope.  At  this  point  the  transducer 
setup  is  complete  and  the  measurement  procedure  is  started. 

The  exact  fundamental  frequency  of  the  crystal  is  first  evaluated 
by  changing  the  frequency  on  the  VFO  until  minimum  attenuation  of  the 
last  echoes  is  observed.  The  frequency  is  then  read  on  the  frequency 
counter.  The  length  of  the  pulse  is  Increased  until  overlap  of  success¬ 
ive  echoes  occurs.  The  frequency  of  the  VFO  is  then  changed  by  about 
1  MHz  on  either  side  of  the  crystal  fundamental  frequency  (fQ).  As  the 
frequency  changes.  Interference  minima  and  maxima  are  observed.  The 
gate  of  the  boxcar  integrator  is  positioned  at  the  overlap  of  the 
second  and  the  third  echo.  The  number  of  minima  between  the  highest  and 
the  lowest  frequency  is  then  noted  and  the  values  of  these  two  fre¬ 
quencies  are  read  from  the  frequency  counter.  The  velocity  v  is  given 
by: 

v  *  2  x  L  x  4f 


Vor  the  velocity  measurements  the  RF  was  set  at  4  MHz  which  was 
approximately  the  fundamental  frequency  of  the  crystal  used  for 
measurements . 


so 


where  L  is  the  sample  length  and  Af  is  the  change  In  frequency  between 
two  "in  phase"  conditions  of  the  second  and  third  echo.  To  improve  the 
accuracy  of  the  velocity  value  21  sets  of  measurements  of  about  15 
minima  each  are  taken.  The  average  velocity  (v)  obtained  is  0.5840  cm/ 
nsec  with  a  standard  deviation  of  0.53».  An  approximate  calculation  of 
the  effect  of  the  phase  change  caused  by  the  relative  magnitude  of  the 
crystal  and  the  sample  mechanical  impedance  is  shown  in  Appendix  C. 

The  result  is  that  this  correction  in  this  case  is  negligible. 

2.  Attenuation  Measurements 

The  procedure  used  for  attenuation  measurements  is  the  same  as 
the  one  used  for  velocity  measurement  up  to  and  including  the  fundamental 
frequency  determination.  The  range  of  usable  frequency  is  3  to  9  MHz. 

The  frequencies  used  are  all  higher  than  3  MHz  to  eliminate  two  problems. 
First,  an  impedance  mismatch  between  the  pulser  and  the  transducer 
(this  mismatch  is  an  inverse  function  of  the  frequency);  and  second,  an 
increasingly  difficult  mechanical  coupling  of  the  crystal  to  the  sample, 
as  the  thickness  of  the  crystal  Is  Increased.  No  frequencies  higher 
than  9  MHz  can  be  used  as  the  Increase  of  attenuation  of  the  ultrasonic 
wave  in  the  sample  prevented  measurement  of  a  sufficiently  large  number 
of  echoes. 

Once  the  resonant  frequency  Is  defined,  the  pulse  parameters  are 
adjusted  to  give  pulses  and  echoes  with  minimum  distortion  (flat  top) 
and  which  show  no  sign  of  Interfering  with  each  other.  The  width  of  the 
gate  Is  then  set  so  that  It  could  cover  only  the  flat  top  part  of  any 
echo.  The  first  echo  is  then  set  at  a  given  reference  level 
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(2  or  5  volts)  and  all  the  other  echoes'  amplitude  measured  by  the 
boxcar  Integrator.  The  reference  setting  Is  checked  regularly  to 
ensure  the  validity  of  the  relative  amplitude  measurements.  The  trans¬ 
ducer  and  the  crystal  are  then  removed  and  the  measurement  repeated  on 
another  facet. 


4*T 


4 
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CHAPTER  IV 


RESULTS  AND  CONCLUSIONS 
A.  EXPERIMENTAL  RESULTS 

1 .  The  Amplitude  Pattern 

A  series  of  amplitude  measurements  was  made  in  accordance  with 
the  procedure  described  in  Chapter  III.  These  measurements  were  made 
for  34  different  combinations  of  four  parameters:  (1)  frequency  (f), 
(2)  wedge  angle  (y),  (3)  sample  length  (L),  and  (4)  radius  of  the 
piezoelectric  element  (a). 

Figure  IV-1  shows  an  example  of  exponential  decay  pattern 
observed  on  the  CRO.  These  patterns  were  observed  when  the  wedge  angle 
was  small  enough  and  the  frequency  low  enough.  In  other  cases  in  which 
the  frequency  was  high  enough  and  the  wedge  angle  large  enough  signifi¬ 
cant  fluctuations  about  the  exponential  decay  patterns  were  observed. 

A  photograph  of  one  of  these  fluctuating  amplitude  decay  patterns,  as 
seen  on  a  CRO,  is  shown  in  Figure  IV-2.  A  plot  of  such  fluctuating 
amplitudes  is  shown  in  Figure  IV-3. 

2.  The  Attenuation  Parameter 

The  attenuation  parameter  a  is  the  quantity  that  is  modified  in 
the  course  of  this  study  to  bring  it  closer  to  the  ultimate  goal  a^. 

The  parameter  a  is  given  by  Eq.  (2-6): 

•  ■  17  ’osio'TOT1  • 
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Figure  IV-Z.  Fluctuations  In  the  Echo  Pattern  Observed  at  3  MHz 
=  11.75  x  10-3  Radians. 
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In  our  study  we  will  use  Pg  as  the  amplitude  of  the  first  measured  echo 
and  P(x)  as  the  amplitude  of  the  nth  echo.  Then  Ax,  the  propagation 
path  from  echo  1  to  echo  n,  is  given  by: 

Ax  =  (R(Pn  '  (Roh 

where  (Rg)n  is  the  distance  between  the  center  of  the  receiver  and  the 
source  after  n  reflections  and  is  given  by  Eq.  (2-42).  Values  of  a  are 
calculated  by  the  computer  program  in  Appendix  B  and  are  given  by  the 
variable  DBMP.  When  we  use  the  uncorrected  data,  this  parameter 
represents  the  apparent  attenuation.  Some  plots  of  apparent  attenuation 
versus  n  obtained  at  4  MHz  and  7  MHz  for  various  wedge  angles  are 
shown  in  Figures  IV-4  and  IV-5.  The  plots  show  that  (1)  the  apparent 
attenuation  is  a  strong  function  of  the  wedge  angle  and  the  echo 
number;  (2)  for  higher  frequency  the  attenuation  parameter  becomes 
sensitive  to  smaller  variations  of  wedge  angle  and  echo  number;  and 
(3)  the  echo  trains  exhibit  both  variations  (monotonic  increases)  and 
fluctuations  (somewhat  periodic  changes). 

B.  THE  DIFFRACTION  CORRECTION 

The  diffraction  correction  defined  by  Eq.  (2-21)  can  be  applied 
to  amplitude  values  obtained  on  facets  with  zero  wedge  angles.  A 
corrected  attenuation  value  can  then  be  obtained  by  using  these 
corrected  amplitude  values  in  Eq.  (2-6).  This  corrected  attenuation 
Is  evaluated  by  the  computer  program  in  Appendix  B  and  is  given  by  the 
variable  DBCD.  A  plot  of  both  the  apparent  and  corrected  attenuation 
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Figure  IV-4.  Apparent  Attenuation  at  4  MHz  for  Different 
Wedge  Angles. 
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Figure  IV-5.  Apparent  Attenuation  at  7  MHz  for  Different 
Wedge  Angles. 
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as  a  function  of  the  echo  number  for  data  taken  at  3.34  MHz  is  shown 
in  Figure  IV-6.  A  similar  plot  for  data  taken  at  6  MHz  is  shown  In 
Figure  IV-7.  The  plots  show  that  the  corrected  attenuation  has  a 
smaller  magnitude  than  the  apparent  attenuation.  For  the  first  two 
echoes  at  6  MHz  and  for  the  first  six  echoes  at  3  MHz  variations  of 
apparent  attenuation  are  larger  for  smaller  frequencies  than  for  larger 
frequencies  in  contrast  with  the  effect  of  the  change  in  wedge  angle. 
(The  change  in  wedge  angle  resulted  In  larger  variations  for  higher 
frequencies.)  These  variations  are  corrected  for  by  the  diffraction 
correction  and  a  reference  attenuation  is  plotted  by  taking  the  average 
of  these  first  few  corrected  echoes.  Dotted  lines  giving  the  reference 
attenuation  are  given  in  Figure  IV-6  and  succeeding  figures. 

Deviation  of  the  data  from  the  reference  attenuation  for  the 
later  part  of  the  echo  train  cannot  be  explained  by  this  correction 
model  for  two  reasons.  First,  even  though  the  approximations  of  the 
model  should  be  better  and  better  as  n  increases  (Rogers  and  Van  Buren 
(1974)  report  an  accuracy  of  at  least  0.6%  in  the  Fresnel  zone  for  the 
range  of  values  of  ka  used  here),  the  deviation  of  the  corrected 
attenuation  from  the  reference  attenuation  becomes  larger  and  larger 
as  n  increases.  Second,  in  contrast  with  the  prediction  of  this 
correction  model,  ths  deviation  of  the  attenuation  data  from  the 
reference  attenuation  is  larger  for  high  frequency  than  the  deviation 
of  the  data  observed  for  low  frequency.  The  possibility  of  having 
Inhomogeneity  In  the  sample,  which  would  give  a  pattern  similar  to  non- 
parallelism,  was  reported  by  Truell  et  al.  (1969).  However,  we  believe 
that  the  quality  of  the  sample  is  such  that  this  possibility  seems 


o(dB/cm) 


a(dB/cm) 


Figure  IV-7.  Attenuation  at  6.00  MHz  for  0  Wedge  Angle. 

A  Reference  Attenuation  Is  Obtained  by  Averaging  the  First  Two  Echoes. 
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unlikely.  However,  this  pattern  Is  very  similar  to  the  results  of 

varying  wedge  angles.  The  accuracy  of  wedge  angle  measurements  was 
-5 

limited  to  3  x  10  radians.  If  we  now  assume  that  a  wedge  angle  of 
2  x  IQ'5  radians  Is  present  here,  then  (from  Eq.  (2-22))  after  three 
reflections  the  Incident  angle  of  the  wave  on  the  receiver  will  have  a 
magnitude  comparable  to  the  smallest  wedge  angle  In  the  sample.  A  con¬ 
clusion  which  could  then  be  drawn  from  the  behavior  of  the  later  part 
of  the  echo  train  is  that,  given  enough  reflections,  the  data  seem  to 
be  more  sensitive  to  wedge  angle  variations  than  the  optical  measure¬ 
ment  method  described  In  Chapter  II. 

These  studies,  then,  allow  us  to  make  two  important  statements. 
First,  for  extremely  small  angles  the  early  part  of  the  echo  train 
allows  one  to  conclude  that  the  diffraction  correction  is  valid  for  up 
to  six  echoes  at  frequencies  as  low  as  3  MHz,  and  up  to  two  echoes  for 
frequencies  as  high  as  6  MHz.  Second,  the  later  part  of  the  echo  train 
follows  the  behavior  predicted  for  any  wedged  sample:  vn.  the  effective 
Incident  angle  of  the  receiver.  Increases  at  each  reflection  so  that  for 
higher  frequencies  the  smaller  wavelength  makes  the  deviation  from  an 
exponential  echo  train  more  noticeable. 

C.  APPLICATION  OF  THE  PUNE  WAVE  GEOMETRICAL  CORRECTION 

It  Is  clear  from  Figures  IV-4  and  IV-5  that  the  magnitude  of  the 
apparent  attenuation  is  strongly  influenced  by  the  value  of  wedge  angle. 
Since  the  diffraction  correction  does  not  take  Into  account  the  change 
in  wedge  angle,  an  attempt  Is  made  to  solve  this  problem  by  the  use  of 
the  Truell  and  Oates  correction  function  H(yn)  given  in  Eq.  (2-27). 
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The  amplitude  values  are  first  corrected  for  diffraction.  From  these 
values,  a  diffraction  corrected  attenuation  (given  by  DBCD  In 
Appendix  B)  is  obtained.  T*en  the  diffraction  corrected  amplitude 
values  are  recorrected  by  using  the  Truell  and  Oates  correction 
function.  The  attenuation  obtained  from  these  recorrected  values  (given 
by  DBCPG  in  Appendix  B)  will  be  referred  to  as  attenuation  corrected 
for  plane  wave  geometry.  The  plane  wave  correction  model  is  now  applied 
to  fluctuating  patterns  such  as  the  one  shown  in  Figure  IV-3  (p.  55). 

The  diffraction  corrected  attenuation  as  well  as  the  attenuation 
corrected  for  plane  wave  geometry  are  shown  in  Figure  IV-8.  As  can  be 
seen  from  these  plots  when  the  fluctuations  of  attenuation  values  are 
fairly  large,  as  long  as  the  propagation  path  is  less  than  2  Zf,  the 
plane  wave  model  corrects  in  the  right  direction.  However,  in  these 
cases  it  significantly  overcompensates  for  the  fluctuations.  Examination 
of  the  data  reveals  that  this  model  Improves  the  prediction  of  the 
shape  of  the  fluctuations  as  we  get  closer  to  the  Fresnel  zone  or, 
equivalently,  as  the  pressure  distribution  on  the  receiver  approaches 
the  plane  wave  model.  This  correction  was  applied  to  different  combi¬ 
nations  of  frequency  and  wedge  angle  and  In  most  cases  It  overcompensated 
for  large  wedge  angles  and  undercompensated  for  small  wedge  angles. 

Calder  showed  that  his  correction  was  quantitatively 
inappropriate  for  his  measurements  In  solid  helium.  The  results  pre¬ 
sented  above  agree  with  Calder.  Even  though  the  Truell  and  Oates  model 
gives  the  approximate  shape  of  the  observed  fluctuations,  it  is  not 
possible  to  use  it  for  quantitative  corrections  in  our  case. 
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Figure  IV-8.  Attenuation  at  4  MHz  for  a  Wedge  Angle  of 
11.75  x  10-3  Radians  Using  the  True!!  and  Oates  Correction. 


D.  APPLICATION  OF  THE  DIFFRACTED  WAVE 
GEOMETRICAL  CORRECTION 

The  diffracted  wave  geometrical  correction  Is  now  applied  to  the 
uncorrected  farfleld  data  and  the  values  of  a(n)  (given  by  the  variable 
DBCDG  In  Append-*  B)  versus  n  are  plotted.  For  comparison  In  the  far¬ 
fleld  and  to  allow  for  partial  correction  In  the  Fresnel  zone,  the 
diffraction  correction  also  is  given. 

Such  corrections  to  the  data  shown  in  Figure  IV-3  (p.  55)  are 
plotted  In  Figure  IV-9.  The  diffracted  wave  geometrical  correction 
brings  the  values  much  closer  to  the  values  observed  for  smaller 
wedge  angles  In  Figure  IV-4  (p.  57).  In  Figure  IV-9  one  observes  that 
the  correction  gives  more  consistent  results  than  the  two  previous 
corrections  and  brings  the  corrected  values  closer  to  the  values 
obtained  for  smaller  wedge  angles  in  the  region  In  which  fluctuations 
occur.  For  some  specific  combinations  of  relatively  large  frequency 
and  wedge  angles,  such  as  the  one  shown  In  Figure  IV-10  for  data  taken 
at  5.75  MHz,  this  geometrical  correction  successfully  corrects  the 
apparent  attenuation  such  that  the  corrected  attenuation  closely 
approximates  the  reference  attenuation  value  obtained  for  zero  wedge 
angle  at  approximately  6  MHz.  In  general,  for  the  frequencies  In 

.3 

the  4  to  6  MHz  range  and  wedge  angles  smaller  than  3.82  x  10  radians, 
the  correction  brought  at  least  the  first  one  or  two  echoes  of  the  far¬ 
fleld  Into  the  range  predicted  from  the  reference  attenuation  obtained 
for  smaller  wedge  angles. 

In  Figure  IV-11  a  plot  of  the  corrected  data  at  3.13  MHz  for  a 
wedge  angle  of  3.82  x  10“3  are  shown.  As  can  be  seen  from  the  plot. 
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Figure  IV-9.  Attenuation  at  4  MHz  for  a  Wedge  Angle  of 
11.75  x  10“ 3  Radians  Using  Data  Corrected  /or  Diffracted' Wave  Geometry 
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3.82  x  10*3 


re  IV-10.  Attenuation  at  5.75  MHz  for  a  Wedge  Angle  of 
Radians  Using  Data  Corrected  for  Diffracted  Wave  Geometry 


and  Reference  Attenuation  for  y  *  0  and  f  *  6  MHz. 
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Figure  IV-11.  Attenuation  at  3.13  MHz  for  a  Wedge  Angle  of 
3.82  x  10’ 3  Radians  Using  Data  Corrected  for  Diffracted  Wave  Geometry 
and  Reference  Attenuation  for  y  *  0  and  f  *  3.34  MHz. 
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the  corrected  reference  attenuation,  at  least  for  the  first  few  echoes, 
gives  a  close  approximation  of  the  reference  attenuation  line.  The 
correction  also  was  applied  at  4  MHz  to  angles  of  3.82  x  10"3  radians 
or  smaller  and  in  all  cases  the  corrected  values  of  the  first  one  to 
three  echoes  were  within  5X  of  the  reference  attenuation  obtained  at 
lower  wedge  angles. 

In  summary,  the  important  points  are  that  (1)  for  frequencies 
nigher  than  4  MHz  this  model  corrects  partially  for  the  deviation  from 
the  reference  attenuation  when  the  type  of  fluctuations  shown  in 
Figure  IV-9  do  not  occur,  i.e.,  for  relatively  low  wedge  angles;  and 
(2)  if  the  frequency  is  4  MHz  or  lower,  then  our  model  successfully 
corrects  for  geometrical  attenuation  (at  least  for  the  first  few 
echoes)  even  for  wedge  angles  as  large  as  3.82  x  10"3  radians. 

E.  CONCLUSIONS 

From  the  analysis  of  the  data  reported  here  two  main  conclusions 
seem  to  emerge. 

1.  The  large  fluctuations  and  variations  observed  when  either 
the  wedge  angles  or  the  echo  number  were  large  enough  or  the  frequency 
high  enough,  cannot  be  explained  entirely  on  the  basis  of  the  diffraction 
model  used  here.  Even  though  the  neglect  of  the  beam  displacement 
and  of  the  actual  pressure  distribution  seems  to  invalidate  the 
Truell  and  Oates  plane  wave  model  In  our  situation,  the  results  show 
that  the  Truell  and  Oates  function  H(yn)  does,  in  contrast  with  the 
other  corrections,  take  into  account  the  form  of  the  large  variations 
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and  fluctuations  at  least  when  the  shape  of  the  wavefront  approaches  the 
plane  wave  approximation. 

Kinsler  et  al.  (1982)  gave  a  description  of  similar  fluctuations 
and  variations.  They  attribute  these  fluctuations  and  variations  to 
a  receiver  directivity  factor  Hr(e.),  where  ei  is  the  incident  angle 
of  the  wave  on  the  transducer.  According  to  Kinsler,  the  receiver 
directivity  factor  H  (e^)  is  given  by  the  magnitude  of  the 
directivity  function  H(e)  of  the  receiver  when  used  as  a  transmitter. 

A  mathematical  justification  of  this  statement  is  shown  In  Appendix  D. 
Then  if  a  plane  wave  is  Impinging  on  a  receiver  Hr(ei)  will  be  given  by 
the  magnitude  of  the  directivity  function  of  a  piston  source.  If  we  now 
take  the  magnitude  of  the  directivity  function  of  a  piston  source  given 
in  Eq.  (2-20)  we  obtain  the  Truell  and  Oates  function  H(yn)  given  in 
Eq.  (2-26)  since  yn  gives  the  magnitude  of  e^.  Then  in  our  case  the 
value  of  the  Truell  and  Oates  function  seems  to  reside  in  the  fact  that 
it  has  exactly  the  same  form  as  the  receiver  directivity  function  for  a 
plane  wave  impinging  on  a  circular  receiver. 

We  are  then  led  to  the  conclusion  that  when  the  wedge  angle  or  the 
echo  number  is  large  enough  or  the  frequency  high  enough,  both  the 
geometrical  attenuation  and  the  receiver  directivity  significantly 
affect  the  measured  attenuation.  And,  in  our  case,  the  model  of  the 
receiver  directivity  factor  based  on  the  plane  wave  assumption  can  only 
be  used  in  a  qualitative  manner.  The  important  point,  however,  is  that 
the  model  of  the  receiver  directivity  factor  gives  a  plausible  physical 
explanation  of  the  large  fluctuations  and  monotonic  variations  observed 
in  the  echo  train. 


2.  From  the  analysis  of  our  data  it  appears  that  a  geometrical 
correction  should  be  applied  to  all  the  amplitude  measurements  made  on 


wedged  samples  when  the  intrinsic  attenuation  is  desired.  Only 
the  cases  when  the  wedge  angle  is  so  small  that  it  cannot  be  detected  by 
a  method  of  accuracy  comparable  to  the  one  used  here,  is  It  sufficient, 
at  least  for  frequencies  lower  than  6  Wz  and  for  the  first  few  echoes, 
to  use  the  diffraction  correction  for  parallel  faces  as  geometrical 
correction.  In  all  the  other  cases  the  diffracted  wave  geometrical 
correction  should  be  used. 

In  some  cases  the  geometrical  correction  is  the  most  significant 
correction  and  the  receiver  directivity  factor  can  be  neglected.  These 
are  characterized  by  the  fact  that  the  values  of  amplitude  corrected 
for  geometrical  attenuation  gives  an  approximately  constant  intrinsic 
attenuation.  In  our  experiments  they  occurred  mainly  at  frequencies 
lower  than  4  MHz  and  for  wedge  angles  smaller  than  3.82  x  10  radians. 

F.  SUGGESTIONS  FOR  FURTHER  WORK 
1 .  Receiver  Directivity  Factor 

Probably  the  most  significant  improvement  to  the  correction  for 
nonparallelism  would  be  the  definition  of  a  receiver  directivity  factor 
which  could  be  applied  to  cases  for  which  the  plane  wave  assumption  is 
not  satisfied.  As  a  secondary  investigation,  a  calculation  of  the 
directivity  function  in  the  farfield  of  a  circular  source  having  a 
distribution  of  the  form  2J^(r)/r  was  made.  The  results  show  that  a 
considerable  decrease  in  the  amount  of  overcompensation  for  the  receiver 
effects  is  obtained  when  this  directivity  function  is  used  for 


correction.  Further  investigation  along  this  direction  probably  would 
be  f rui tf ul . 


2.  Correction  in  the  Fresnel  Zone 

A  very  useful  improvement  of  the  diffracted  wave  geometrical 
correction  would  be  to  extend  its  validity  farther  into  the  Fresnel  zone. 
This  would  allow  for  correction  in  a  region  which  is  commonly 
encountered  in  nondestructive  testing  situations.  The  mathematical 
difficulty  associated  with  this  improvement  is  a  long-standing  impediment 
to  solution  of  all  diffraction  problems. 

3.  The  Gaussian  Transducer  Model 

During  the  course  of  this  investigation  it  became  clear  that  an 
analytical  solution  based  on  a  radiating  element  with  a  diffraction 
field  showing  radial  symmetry  would  be  difficult  to  handle.  This  is  the 
reason  a  numerical  solution  was  used.  An  experimental  configuration 
which  would  lend  itself  to  the  use  of  cartesian  coordinates  would  be 
much  simpler  to  model  mathematically.  The  Gaussian  transducer  developed 
by  Martin  and  Breazeale  (1971)  has  this  characteristic. 

A  theoretical  solution  of  the  problem  of  evaluating  the  diffraction 
field  of  a  Gaussian  transducer  is  included  in  Appendix  E.  Experiments 
modeled  on  this  solution  would  be  very  interesting. 

4.  Velocity  Measurements 

As  can  be  seen  from  Figures  I 1-9  and  11-10  (pp.  26  and  29),  if 
the  wedge  angle  or  the  sample  length  is  large  enough  significant 
variations  in  the  measurements  of  velocity  will  occur  because  of  the 
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variations  in  the  propagation  path  between  echoes.  Preliminary  results 
obtained  in  large  steel  blocks  (20  cm)  with  large  wedge  angles  (1  to  2°) 
at  a  frequency  of  about  5  MHz  showed  that  significant  decreases  (up  to 
5%)  can  be  observed  in  the  transit  time  between  two  successive  echoes. 

A  study  of  the  variation  of  velocity  measurements  in  wedged  samples 
promises  to  be  an  interesting  field  of  investigation. 
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APPENDICES 


APPENDIX  A 


AN  APPLICATION  OF  THE  RECIPROCITY  THEOREM 
A.  THE  RECIPROCITY  THEOREM 

Let  us  assume  two  identical  transducers  1  and  2.  Let  a  volume  V 
with  closed  surface  S  bound  these  two  sources  but  not  enclose  them,  as 
shown  in  Figure  A-l.  ^  is  the  velocity  potential  for  transducer  1 
and  <|>2  for  transducer  2.  If  we  excite  transducer  1  only,  the  description 
of  the  oropagation  of  the  wave  in  volume  V  is  given  by: 

*  -K241  (A-l) 

and  the  particle  velocity  u1  =  v<|>^  and  the  pressure  P-j  *  -ju>p0<|>.|. 
Similarly,  if  we  now  excite  transducer  2  only  at  the  same  frequency  u>, 
we  have: 

v2*2  =  -K202  (A-2) 

and  again,  u2  ■  v<t>2  and  P2  *  jup04>2. 

From  Green's  theorem: 

fs  (*ly*2  “  Vh*  *  "  dS  *  |v  “  V^l  )dV  •  (A'3) 

Substituting  (A-l)  and  (A-2)  in  (A-3)  we  have: 

-jwpQ  f  (P,  U2  -  P2  u})  •  n  dS  -  |  (-^K2^  ♦  ♦jA^dV  .  (A-4) 
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The  right  hand  side  of  (A-5)  vanishes  identically  and  we  have: 


n  dS  *  0  . 


This  is  the  expression  of  the  principle  of  reciprocity.  Let  us  now 
restrict  this  theorem  to  the  situation  where  the  closed  surface  is  made 
of  an  infinite  rigid  baffle  of  surface  Sg  such  that  (u*n)  *  0  on  Sg. 
Then  the  only  surface  allowed  to  react  to  the  pressure  waves  are  the 
two  transducer  surfaces  and  Sg.  Then,  (A-5)  becomes: 


I  <pi  “2  -  P2  fy  *  n  dS  =  0  (A-6) 

^ (s1+s2) 

or  equivalently. 


(A-7) 

Since  the  original  assumption  is  that  one  transducer  is  passive  while 
the  other  is  active,  then  wh^n  one  transducer  is  active  the  pressure 
over  the  surface  of  the  second  transducer  is  not  zero  but  the  velocity 
over  that  surface  is  zero.  This  means  that  the  negative  terms  of 
Eq.  (A-7)  go  to  zero  giving: 


(A-8) 
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B.  SMALL  SOURCES  IDENTITY 


If  the  two  sources  1  and  2  are  small  with  respect  to  the 
wavelength  and  located  several  wavelengths  apart  such  that  the  pressure 
is  uniform  over  each  source,  then  (A-8)  when  time  is  considered  can  be 
expressed  as: 


jr]t7  J  u,(t)  •  "  <ts  •  frjtT  f  :2(t)  •  " 


(A-9) 


We  now  define  the  rate  at  which  the  volume,  in  the  vicinity  of 


the  source,  is  displaced  in  function  of  the  source  strength  Q  as: 


Q  e^ut  -  j  u(t)  •  n  dS  , 


(A -10) 


Using  the  fact  that  the  variation  of  P  is  given  by  P-j(t)  =  P-j  eJ(1)t  and 
Eq.  (A- 10),  we  obtain  from  Eq.  (A-9)  that 


Ql  Q2 


(A-ll) 


Therefore  the  ratio  of  source  strength  to  pressure  for  a  small  source 
is  independent  of  the  shape  of  the  source. 


n  n  tin  nun  pnnn  luiiinnnnn  nm*nm»n 


APPENDIX  B 


A  COMPUTER  PROGRAM  FOR  CORRECTING  ULTRASONIC 
ATTENUATION  DATA 


THIS  PROGRAM  IS  CODED  IN  FORTRANIO  LANGUAGE. 

IT  CORRECTS  AMPLITUDE  VALUES  OBTAINED  IN  A  WEDGED 
SAMPLE  WITH  A  CIRCULAR  ULTRASONIC  TRANSDUCER. 

IT  ALSO  CALCULATES  THE  CORRECTED  ATTENUATION 
BETWEEN  ECHO  1  AND  N  I N< 151 


THE  INPUT  VARIABLES  ARE: 

VS  VELOCITY  OF  THE  WAVE  IN  THE  SAMPLE  (CM/USFCI 
F:  FREQUENCY  (MHZ» 

D:  DIAMETER  OF  THE  PIEZOELECTRIC  ELEMENT  CCMI 
SL:  SAMPLE  LENGHT  ICMl 
GAMMA:  WEDGE  ANGLE  (RAOIANSi 

THE  OUTPUT  VARIABLES  ARE  : 

AXXX :  AMPLITUDE 

DB XXX:  ATTENUATION  (DEC IBELS/CM) 

WHERE  XXX  MEANS: 

MP:  UNDIRECTED 

CD:  DJ  ‘'.ACTION  CORRECTED 

CDGi  CORRECTED  FOR  DIFFRACTED  WAVE  GEOMETRY 
CPS:  CORRECTED  FOR  PLANE  WAVE  GEOMETRY 


IMPLICIT  DOUBLE  PRECISION  CM) 

DIMENSION  AMPC15I  ,ACD(15),ACPG(15  ),ACDG(  15) 
DIMENSION  RMPCI5)  »RCD(  15)  »RCPG(  15  ),RCDG(  15) 
DIMENSION  RCC15  )»ZC(  151 

COMMON  P I  *  4  *  SL ,GAM2N»WK ,OME  G A, ALAMBD, ECHONO 
COMMON  ROD , ZO, S, SRF 
PI  »  3.1415  926 
V  «  5.840E-I 
10  READ  (1,120)  F,  D,  SL, GAMMA 

A  »  D/2. 

IF  (F.GT.100.01  STOP 

ALAMBD  «  V/F 
OMEGA  -  2.*PI*F 
MK  *  2**PI/  ALAMBD 
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WRITE  (5,130)  A, SI, GAMMA, F 
ICOUNT  *  1 

READ  (1,10)  AMPt  ICOUNT  ) 

IF  (AMP( ICOUNT). ST. 100.0)  GO  TO  30 

ICOUNTMCOJNT+1 

GO  TO  20 

N  *  ICOUNT- 1 

DO  90  I  *  1  ,N 

ECH0N3  *  0s  LOAT  ( I  ) 

GAM2N  *  2.0*GAMMA*ECH3N3 
IP AR  *  ECH0NO*$L*2.0 
S*  ( ZP AR* ALAMBD) /( ( A )**2 ) 

CALL  DIFCOR  (CO) 

ROO  *  0.0 
20  «  2PAR 

AC D(  I)*(AMP(I) )/(CD) 

CALL  FLUCTJ(HGAM2N) 

FARFLO  *  (0**2-ALAMBD**2)/(4*ALAMBD) 

IFIZPAR.LT. { 1 . 0*F ARFLO) )  GO  TO  60 

CALL  DSCOR  I  COG) 

ACDGU)  *  AMP(  I  )/CDG 
RC ( I )  «  R3D 
ZC(I  )  *  ZD 

ZF  «  SQRT(ROO**2*ZO**2) /FARFLO 
IF  (ZF.LT.3.0)  GO  TO  70 

GO  TO  80 
ACDG( I )  *  A  CO (  I  ) 

ACPG(I)  *  ACOG( Il/HGAWPN 
RC  ( I  )  =0. 

ZC(I)  *  ZPAR 

GO  TO  90 
CONTINUE 

WE  NOW  CORRECT  THT  FARF IELO  APPROXIMATION  TO  IMPR 
ITS  CONVERGENCE • 

IF (ZF.LT .1. 5)CF*( ZF- 1. ) *0.6*0. 4 

IF  (ZF.GT.1.5.AND.ZF.LT.2.)CF«(ZF-1.5)*0.4*0.7 

IF  (ZF.GT. 2.0.  AND.  ZF.LT  • 2. 5) CF* (  ZF— 2  •  )<  0.14*0.9 

IF(ZF.GT.2.5)CF«(ZF-2.5)*0.04*0.97 

ACDGU)  *  ACDG(  I)  *CF 

ACPG(I)  «  ACD( I J/HGAM2N 

CONTINUE 

WRITE  (5,150) 

DO  10D  I  «  1 , N 

RMP( I )  «  AMP( I ) /AMP ( 1) 

RC  0( I )■  ACD (I  )  / ACO( 1) 

RCDGU  )  -  ACOG(  l)/ACDG(  1 ) 

RC PG ( I ) *  A:PG(I)/ACPG( l) 

X*SQRT(KC(! )**2  *ZC( I )**2 )-SQRTI RC 1 1 ) **2*ZC I 1)**2) 
WRITE  (5,180)  I,AMP(I),ACD(1),AC0G(I),ACPG(I) 

IF  (  I.EQ.l)  GO  TO  100 


IMPROVE 
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DBMP  *  2 0. 0*ALOG10 l RMP ( I ) )/X 
06 CD  ■  23 .0*AL0G1 0(RCD( I) I  /X 

OBCOG  >  20. 0* A LOG 10 (RCDG(I))  /X 
DBCPG  «  20.0*AL0G10(RCPG(I) I  /X 
WRITE  (5*170)  DBMP, DBCD* OBCOG >  OBCPG 
100  CONTINUE 

GO  TO  10 
RETURN 

110  STOP 

120  FORMAT  { 3F8 , 8X, El  1. 3 ) 

130  FORMAT  (1H1 »4(1X»1PE13.5)I 

143  FORMAT  (  F5 I 

150  FORMAT  (IX, 'THE  RESULTS  ARE  GIVEN  AS '/IX, ‘ECHO  NUMBER 

1  AMP'tlSX,*  ACD ' , 1 6X  * ' ACOG' ,  16X  ,* ACPG* /12X,* 

2  DBMP* , 15X* • DBCD* ,15X, ‘OBCOG I5X, • DBCPG* I 

160  FORMAT  (lN-.'ECHO  1  *  *  IX ,4 ( 1PE1 6. 9) ) 

170  FORMAT  (  1H0 , 8X  ,  4(  1PE 16.  9H 

180  FORMAT  ( IH~  »' ECHO  • *  12, 1 X, 4 ( 1  PE  16 .9 ) ) 

END 

C 


SUBROUTINE  OIFCOR(CO) 

THIS  IS  ROGERS  AND  VAN  BUREN  CORRECTION  FOR  DIFFRACTION 


IMPLICIT  DOUBLE  PRECISION  (Ml 
COMMON  PI,A,SL,GAM2N,WK, OMEGA, ALAMBD.ECHONO 
COMMON  R03  ,ZO,S,SRF 
MARG«2.*PI/S 

G  CALL  THE  IMSL  FUNCTIONS 

DJO  *  MMBSJ  OCMARGtIER) 

DJI  *  MMBSJ  KMARG,  IERI 

10  CD  *  S3RT( ( 0  OS ( MARG) -DJO )**2  +  (SINC  MARG) -DJI )**2I 

RETURN 
END 


SUBROUTINE  FLUCTU  ( HGAM2NI 

THIS  IS  TR J ELL  AND  OATES  FLUCTUATION  CORRECTION 


IMPLICIT  DOUBLE  PRECISION  (M) 

COMMON  PI , A , SL, GAM2N,WK* OMEGA, ALAMBD.ECHONO 

COMMON  ROO ,Z0, S, SRF 

MARG  *  ABStWK*A*SINI  GAM2N) ) 

HGAM2N  «  1.  0 

IF  (  MARG.LT .1.00-8)  GO  TO  10 


L>  UIU.HI  DO 
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:  CALL  THE  IMSL  FUNCTION  MMBSJi 

OJI  «  MMBSJ 1(MARG,IER) 

HGAM2N  «  ABS(2.0*DJl/MARGI 
10  CONTINUE 

RETURN 
END 

m 
< 0 

C 

SUBROUTINE  DSKCOR  (CDGI 

Z  THIS  IS  THE  DIFFRACTED  WAVE  GEOMETRY  CORRECTION 


IMPLICIT  DOUBLE  PRECISION  (M) 

COMMON  PI ,A,SL,GAM2N,WK, OMEGA, ALAMBD.ECHONO 

COMMON  RDO  »Z0, S* SRF 

EXTERNAL  F ARP , FARPD 

ROO  *  ECH0N0*SL*SIN(GAM2N) 

TAN  *  S INI  3 AM2N/2 •  I /CDS  I GAM2N/2 . ) 

zo  *  2.*e:hono*sl-roo*tan 

CALL  THE  IMSL  NUMERICAL  INTEGRATION  FUNCTION  DBLINT 
WHERE  FARP  AND  FARPD  DEFINE  THF  PRESSURE  AT  A 
POINT  (RtTHETA)  LOCATEO  IN  THE  FARFIELD. 


R1  *  0.0 
R?  *  A 

AERR  *  l.OE-3 
THETA1  *  0.0 
THETA2  *  2 •  *PI 

C* DBLINT  (FARP,Rl,R2,THETAl ,THE TA2 , AERR, ERROR , I ER I 
CD=DQL INT(FARPD»R1,R2»THETAI» THET A2, AERR, ERROR, IERI 
SRF  *  PI  *(  A **2) 

COG  =  SORT! (C/SRF J**2*l CD/SRF 1**21 

RETURN 

END 


REAL  FUNCTION  FARP  IR.THETAI 
IMPLICIT  DOUBLE  PRECISION  IN) 

COMMON  PI,A,$L, GAM2N,WK , OMEGA, AL AMBD. ECHDND 
COMMON  ROOtZO, S, SRF 

RO  «  ABS ( SORT (( (ROO  +R*COS (T HETA l*COS (GAM2N ) 
II **2)*(RPSIN(  THETA  11**211 
Z  «  ZO  -  R*  COS ( THET  A )*S INIGAM  2N ) 

DIR  *  1. 

OROZ  ■  SORT  (R0**2*Z**2) 

MARG  «  WK*A*RO/ORQZ 
OJI  >  MMBSJ  i(  MARG ,  I  ER  I 
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IFIMARG.GT. t .0-8) 0IR*2- *DJ1/NARG 

PA  XI  $*2*SIN(  <  PI  /ALAMBO)  *  (  SQRTI  A**2«-0R0Z**2)-3RDZ  1 1 
FARP  ■  OIR*COSIW<*QROZI*R*PAXIS 
10  CONTINUE 

RETURN 
END 

REAL  FUNCTION  FARPD  IR « THETAI 
IMPLICIT  DOUBLE  PRECISION  I M) 

COMMON  PI ,A,SL,GAM2N,WK, OMEGA, ALAMBD.ECHOND 
COMMON  ROO  »  Z0»  S, SRF 

RO  »  ABS ( SQRT ( ( (ROO ♦R*CO S ( THETA )*COS (GAM2N) 

11  **2 1  *<R*S  IN  I  THETA  H**2>) 

Z  *  ZO  -  R*COS( THET A )*S INI 3AM2N I 

OROZ  «  SORT (RO**2+Z**2> 

01 RO  =  1. 

MARG  *  WKOA*R0/0R0Z 

OJl  »  MMB SJ 1 < MARS  « I ER) 

IFIMARG.GT. 1 • 0-8  1 01 RD*2.*0J l/MARG 

PAXI S=2*SINI (PI/ALAMBDI *1 SQRT I A**2+OROZ**2 1-3R3Z 1 1 
FARPO  »  0IR0*SIN( WK*OROZI*R*PAXIS 
10  CONTINUE 

RETJRN 
END 


APPENDIX  C 


PHASE  CORRECTION  FOR  VELOCITY  MEASUREMENTS 

As  pointed  out  by  many  authors  (McSkimin,  1950;  Williams,  1958; 
Papadakls,  1967;  Truell,  1969),  when  a  wave  Is  reflected  from  an  inter¬ 
face  a  change  of  phase  occurs.  The  phase  shift  is  a  function  of  the 
relative  velocity,  and  the  mechanical  impedance  of  the  two  materials  in 
contact  at  the  interface  as  well  as  the  frequency  of  the  wave.  This 
means  that  we  have  to  correct  for  two  Interfaces:  the  sample- 
couplant  interface  and  the  coupl ant-crystal  interface.  However,  as 
mentioned  by  Williams  (1958)  and  Papadakis  (1967),  the  effect  of  the 
bond  material  can  be  neglected  when  quartz  crystals  of  frequency  of 
10  MHz  or  lower  are  wrung  on  the  surface  of  a  solid.  Therefore,  the 
only  Interface  of  concern  here  is  the  crystal-sample  interface.  A 
correction  AT  for  the  measured  transit  time  T  was  given  by  Papadakis 
(1967)  for  a  single  transducer  pulse  echo  overlap  method  which  is  very 
similar  to  the  method  used  here.  Neglecting  the  coupling  layer,  the 
condition  for  overlap  gives 

T  -  l/fQ 

where  fg  Is  the  resonance  frequency  and  the  correction  for  transit  time 
in  the  bond  is 


AT|  «  Yi/2irfj 
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where  y^  Is  the  phase  shift  given  in  radians  for  frequency  f-j .  As 
pointed  out  by  Williams  (1958),  the  factor  y  changes  signs  depending  on 
whether  the  frequency  f  is  higher  or  lower  than  the  resonant  frequency 
fp.  Furthermore,  when  the  coupling  layer  Is  negligible,  y  is 
proportional  to  the  frequency.  Therefore,  by  dividing  the  correction 
AT  into  two  parts,  one  for  frequencies  higher  than  fp  (aTq)  and  one  for 
frequencies  lower  than  fp  (aT^)  and  taking  the  average  (aT)  of  both 
(since  about  the  same  number  of  minima  were  observed  on  both  sides  of 
fp) ,  a  good  indication  of  the  magnitude  of  AT  can  be  obtained.  That  is: 

AT  =  (ATL  +  ATH)/2 

where 

ATL  =  +Y|_/2irfL 
ATh  *  +YH/27rfH 

and 

yh  =  +|2  tan-1  a | 
yl  *  -|2  tan"1  a] 


zt 

a  *  p-  [tan  et]  . 


where 
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The  factor  z  refers  to  the  mechanical  Impedance  and  the  Indices 
t  and  s  refer  respectively  to  the  transducer  piezoelectric  material1 
and  the  sample.2  The  factor  et  is  given  by: 

6t  *  2  X  n  X  f  X 

where  a  is  the  crystal  thickness  and  v  the  velocity  of  propagation  of 
the  ultrasonic  wave  In  the  crystal. 

Applying  the  correction.  It  was  found  that  the  time  correction 
aTl  applied  for  the  average  low  frequency  fL  of  the  21  sets  of 
measurements  was  approximately  -0.1%,  whereas  the  correction  of  the 
time  due  to  the  frequencies  higher  than  the  center  frequency  aTh,  using 
the  average  high  frequency  fH,  was  approximately  +0.12%.  The  result 
of  these  two  corrections  Is  that  they  almost  exactly  compensated  for 
each  other  and  the  average  correction  AT  is  about  50  times  smaller  than 
the  standard  deviation  of  the  measurement. 

The  values  of  the  parameters  used  are  listed  in  Table  C-l. 


Vhe  values  of  z*,  *t»  and  vt  were  obtained  from  Krautkramer 

(1969). 

2zs  was  obtained  by  evaluating  the  density  ps  of  the  sample  and 
using  the  equation  zs  ■  vs  x  ps  where  the  previously  uncorrected 
velocity  v  was  used  as  vs. 
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TABLE  C-l 

PHASE  CORRECTION  DATA 


fH  =  4.550  MHz 
fL  =  3.580  MHz 
f0  *  4.003  MHz 
zt/zs  *  0.3200 

=  0.07175  cm 
*  0.56984  cm/wsec 
yl  =  -.2040 
y^  =  +.3038 
T  =  8.5042  psec 
aTl  *  -9.069  x  10" 3  psec 
A^  *  +10.626  x  10  3  psec 
AT  *  7.783  x  10’4  psec 


AT/T  «  0.009% 


APPENDIX  D 

THE  RECEIVER  DIRECTIVITY  MODEL 

The  magnitude  of  the  mechanical  deformation  of  the  piezoelectric 
element  Is  dependent  upon  the  distribution  of  the  pressure  vectors  on  the 
receiver  surface.  If  one  assumes  that  an  Infinite  plane  wave  Is  Impinging 
on  the  receiver,  then  the  magnitude  of  the  pressure  deformation  on  the  j 

receiver  |P|r  can  be  related  to  the  average  maximum  pressure  of  the  wave 
P  by  the  use  of  a  directivity  function  Hr(e^)  where  e-  is  the  incidence  j 

i 

angle: 


|P|r  -  P  Hr(0i)  .  (D— 1 ) 

Klnsler  et  al.  (1982)  report  different  applications  of  the 
directivity  function  and  show  mathematically  that  the  reciprocity  theorem 
can  be  Invoked  to  prove  that  the  directivity  factor  of  a  receiver  is 
given  by  the  magnitude  of  the  directivity  function  of  the  same  vibrating 
element  when  used  as  a  transmitter.  A  mathematical  proof  as  applied  to 
a  plane  wave  propagation  In  shown  here. 

We  now  assume  the  sources  1  and  2,  described  In  Appendix  A,  to  be 
large  compared  to  the  wavelength  and  made  of  an  active  rigid  piston 
located  In  an  Infinite  rigid  baffle  such  that: 

u  •  n  *  u  (D-2) 

then  Eq.  (A-8)  becomes 
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u2  <(P)r  S2  (S2>  ■  U1  <tP)2>  S1  (S1}  (°-3) 

where  <(P)^>  ^  means  the  pressure  P  Integrated  on  the  surface  Sj  when 

J 

only  the  transducer  1  is  active.  Furthermore,  we  will  assume  that:  the 
two  transducer  acoustic  axes  both  make  an  angle  e  with  the  wave  coming 
from  the  other  transducer;  the  distance  R  between  the  center 
of  the  two  transducers  is  large  enough  that  Fraunhofer  diffraction 
conditions  exist;  andthe  ratio  is  so  small  that  the  wave  originating 
on  S.  appears  plane  when  it  reaches  the  transducer  S^.  This  is  shown 
in  Figure  D-l.  The  pressure  felt  on  element  2  due  to  the  source  1 
<(P)l>  s  can  be  defined  as  a  function  of  the  pressure  field  of  P1 
called  <P,(R,e)>c  where  e  is  the  angle  between  the  acoustic  axes. 

Since  we  have  plane  waves,  then  this  pressure  field  is  constant  over 
the  surface  S2.  Then: 

<pi (r.®)>s2  =  pi(0)  x  A(R)  (D'4) 

where  A(R)  is  a  pressure  amplitude  function.  Similarly 

<(p)2>S]  *  P2(e)  x  A(R)  .  (D-5) 

Replacing  (D-4)  and  (D-5)  in  (D-3)  gives 


A(R)  u2  P-j (e)  S2  *  A(R)  u1  P2(e)  . 


(D-6) 
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In  Figure  D-2  both  acoustic  axes  of  the  two  transducers  are  coincident, 
then  (D-6)  transforms  to: 

#(R>  “2  \x1s  *2  ■  «»)  “1  P2axU  A,  .  (0-7 

The  magnitude  of  the  ratio  of  (A-16)  over  (A-17)  gives 


P^e.t) 

P2(e,t) 

P, 

axis 

^2 

^axis 

If  we  now  define  a  directivity  factor  H^(e)  as 


H^e)  = 


Me) 


axis 


(D-9) 


then  replacing  (A-18)  in  (A-19)  we  have 

H-j  (e)  =  H2(e)  .  (D-10) 

Since  this  equation  must  be  true  independent  of  which  transducer  is  a 
transmitter  and  which  is  a  receiver,  and  since  both  transducers  are 
Identical,  then  we  must  conclude  that  the  receiver  directivity  factor 
H(e)  is  the  same  as  the  transmitter  directivity  factor. 


5 

£ 
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APPENDIX  E 


THE  GAUSSIAN  TRANSDUCER 


Let  us  assume  that  a  pressure  disturbance  has  a  Gaussian  form 
described  by  the  relation: 


P(p,0)  =  PQ  e 


2  /t  2 

P  /Tq 


(E-l) 


where  p  is  the  radial  coordinate  of  a  cylindrical  coordinate  system 
originating  at  the  center  of  the  source.  TQ  is  the  radius  for  which  the 
amplitude  decreases  to  1/e  of  its  maximum  value.  Now,  let  this  disturb¬ 
ance  propagate  in  a  lossless  fluid  a  distance  z.  Since  the  Fourier 
transform  of  a  Gaussian  distribution  is  a  Gaussian  distribution  the 
pressure  P(p,z)  should  have  the  form  of  Eq.  (E-l). 


P(p,z)  =  P(z)  e 


Then  from  the  principle  of  conservation  of  energy: 


(E-2) 


r  p(p.o)  *  r  p(p,z)  , 

Jo  jo 


(E-3) 


or  equivalently. 


2  tj  2 

■P  /Tq 


dp  =  r 

Jo 


-p  /T(z) 


(E-4) 
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where  ?z  is  the  pressure  on  the  axis  of  propagation  at  a  distance 
the  origin.  Equivalently  then, 


From  Abramowitz  (1972): 


r 

Jc 


e'^dp 


T  A/2  . 


Then  replacing  (E-6)  in  (E-5) 


From  Haselberg  and  Krautkramer  (1959): 


Pz  -  P0Mz) 

where  *(z)  ■  [(zx/„T02)2  ♦  1]VJ'  The"  rcpl”c,n3  <E-®>  1n  (E-7): 
T(Z)  •  »(Z)  T0  . 


Replacing  (E-9)  and  (E-8)  in  (E-2)  we  obtain: 


P(p’z)  *?^7e 


P«  -P2/*2(z)Tn2 


-ABBfc  •  ■*  r-t-  •  *■*•-**  *■  •'#, »f  .  C1M 


z  from 


(E-5) 


(E-6) 


(E-7) 


(E-8) 


(E-9) 


(E-10) 
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We  now  assume  a  plane  wave  propagation J  Then 


p,  ,  ..  .  p0  -»V(r)T02  1(„t.kl) 
P(p,z,t)  TfzT  e  e 


(E-ll) 


w  the  magnitude  of  the  average  pressure  on  the  receiver  Is  obtained  by 
integration  of  (E-ll): 


-p2/42(z)T02 

e _ _ 

*Tz) 


(E-12) 


Following  the  same  process  described  In  Chapter  II  we  will  apply 
this  unidirectional  wave  propagation  model  to  a  multiply  reflected  wave. 
We  now  assume  that  the  dimension  of  the  receiver  in  the  y  direction  is 
small  enough  that  the  amplitude  does  not  vary  significantly  along  the  y 
dimension  of  the  receiver.  Then  Eq.  (E-12)  becomes: 


-p(x)V(z(x))T02 

dTKTJl 


e-iMx)dx 


(E-13) 


where  Ax  is  the  length  of  the  receiver  electrode.  Now  we  replace  the 
value  of  z(x)  and  p(x)  by  the  coordinate  system  developed  in  Chapter  II. 
Then 


|P(z)| 


n 


PQ  r  AX/2  -(xcos(Yn)-p(0,0)n)2A2(z(0,0)n+xs1n(Yn)) 

4x  J-ax/26 


■  -y— 1 _ _ _ _e-1Mz(0,0)  n+xsin(Yn))rf  ,F  ... 

♦Wfl)n  +  xsin(Yn))e  n  n  dx  * 


^According  to  Haselberg  and  Krautkramer  (1959)  this  introduces  an 

error  on  P  smaller  than  e”a  ^0  where  a  is  the  radius  of  the  piezo¬ 
electric  element. 


The  circular  piezoelectric  element  and  the  mechanical  coupling 


technique  described  In  Chapter  III  are  also  used  with  the  Gaussian 
transducer.  The  main  difference  between  the  Gaussian  transducer  and  the 
previously  described  transducer  Is  that  now  the  source  electrode  Is  a 
rectangular  strip  of  aluminum  wrung  onto  this  piezoelectric  element. 

By  positioning  the  circular  piezoelectric  element  such  that  the 
rectangular  source  electrode  has  its  length  normal  to  the  slope  of  the 
face  It  is  then  possible  to  use  two  smaller  rectangular  electrodes  as 
receiver  by  placing  themon either  side  of  the  source  electrode  such  that 
the  pattern  of  electrodes  looks  like  a  cross.  Then  the  displacement  and 
rotation  of  the  propagating  wave  which  occurs  strictly  in  the  p'z'  plane 
can  be  handled  by  means  of  the  cartesian  coordinates  as  described  above. 
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